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ABSTRACT 


This  work  comprises  a  theoretical  and  experimental  study  of 
macrosegregation  in  solidification  of  Ingots  and  castings.  A 
general  expression  is  given  for  macrosegregation,  which  is  viewed 
as  resulting  from  mass  flow  into,  or  out  of,  local  regions  of  the 
ingot  during  solidification.  The  mass  flow  considered  in  detail 
results  from  flow  of  solute  rich  liquid  to  feed  solidification  and 
thermal  contraction. 

Detailed  solutions  are  given  for  unidirectional  heat  flow  (but 
where  fluid  flow  is  not  necessarily  unidirectional).  It  is  shown 
that  "inverse  segregation"  and  at  least  some  types  of  macroscopic 
"banding"  can  be  quantitatively  understood  in  terms  of  the  analysis. 
Mold-metal  interface  resistance  is  shown  to  strongly  influence 
macrosegregation,  and  variations  in  heat  transfer  rate  during 
solidification  are  shown  to  result  in  high  or  low  solute  bands. 

Changes  in  solidification  cross-section  also  result  in  macrosegregation. 

Solutions  of  the  macrosegregation  equations  are  next  examined  for 
cases  where  heat  flow  is  multidirectional.  It  is  shown  that  center- 
line  segregation  can  result  from  the  same  basic  mechanism  examined 
above;  i.e.,  flow  of  liquid  to  feed  solidification  contractions.  Other 
types  of  segregation  considered  include  (1)  the  "negative  cone"  of 
segregation  often  found  at  the  base  of  ingots  and  (2)  high  solute 
regions  found  under  some  risers  (hot  tops) .  These  types  of  segregation 
also  result  in  whole  or  in  part  from  the  same  mechanism. 

A  subsequent  section  considers  briefly  some  types  of  segregation 
which  may  result  primarily  from  thermal  contraction  or  other  types  of 
mass  flow  of  solid  within  the  solidifying  ingot.  This  mass  flow  of 
solid  can  accentuate  centerline  segregation,  cause  "hot  tearing"  types 
of  segregation,  and  result  in  segregation  from  settling  or  floating 
of  solid  in  early  stages  of  solidification 

Numerical  examples  of  macrosegregation  are  given  for  selected 
limiting  cases;  calculations  are  based  on  simplified  forms  of  the 
general  macrosegregation  equation  to  permit  analytical  solution. 

Results  of  laboratory  experiments  which  qualitatively  and  quantita¬ 
tively  confirm  the  analysis  are  presented.  Results  are  compared 
with  simple  analytical  solutions,  and  with  more  precise  results  of 
numerical  analysis  solutions  based  on  actual  thermal  data. 

The  experiments  reported  are  on  aluminum-4.5  per  cent  copper 
alloy,  and  include  experiments  on  the  following  mold  designs: 

(1)  simple  unidirectional,  (2)  unidirectional  with  restricted  fluid 
flow,  (3)  unidirectional  with  variable  heat  extraction  rate  to 
illustrate  banding,  (4)  unidirectional  with  reduced  cross  sections, 

(5)  edge-risered  unidirectional  to  illustrate  centerline  segregation, 

(6)  ingot  with  bi-directional  heat  flow  to  illustrate  centerline 
segregation. 
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PREFACE 

This  report  summarizes  the  fourth  year  of  research  on  study  of 
effects  of  solidification  variables  on  casting  and  ingot  structure. 

All  experimental  work  has  been  on  aluminum  and  magnesium  base  alloys, 
but  portions  of  the  work  are  equally  applicable  to  ferrous  and  other 
alloys. 

Work  during  the  first  two  years  of  the  research  concentrated 

primarily  on  microsegregation.  Analyses  were  developed  for  prediction 

of  effects  of  alloy  and  solidification  variables  on  microsegregation; 

these  have  shown  excellent  agreement  with  experiment*  ^ .  Work  during 

the  second  year  also  considered  theoretically  and  experimentally  the 

2  5 

problem  of  pore  formation  in  solidification  '  .  Finally,  in  the 

second  year,  consideration  was  also  given  to  dendrite  structure,  and 

2  4  6 

effects  of  solidification  variables  on  this  structure  *  *  . 

In  the  third  year  of  the  research,  detailed  study  was  made  of 
elimination  of  microsegregation  by  solution  treatment^.  This  work  was 
undertaken  because  previous  studies  had  shown  that  mechanical  pro¬ 
perties  of  both  cast  and  wrought  aluminum  alloys  are  strongly  depen- 

8—io 

dent  on  degree  of  homogeneity  obtained  in  the  solution  treatment 
The  previous  studies  also  showed  it  is  extremely  difficult  to  obtain 
a  high  degree  of  homogeneity  in  most  such  alloys.  A  detailed 
understanding  of  the  kinetics  of  solutionizing  aluminum  alloys  was 
therefore  felt  to  be  of  considerable  engineering  importance. 

Note:  References  for  this  Preface  will  be  found  on  the  fol? owing 

References  for  the  rest  of  the  report  will  be  found  at  the  end 
of  the  report. 
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In  the  third  year  of  the  research,  work  was  also  conducted  on  the 
outer  "chill  zone"  found  in  ingot  solidification,  and  it  was  shown  that 
the  classical  theory  of  "copious  nucleation"  m  the  region  is  inade¬ 
quate  to  explain  experimental  results.  Dendrite  structure  at  the 
chill  surface  was  shown  to  be  different  from  that  in  the  bulk  metal 
in  some  rather  surprising  ways^. 

Finally,  a  series  of  experiments  were  conducted  from  which  it  was 
concluded  that  dendrite  arm  spacing  in  usual  casting  and  ingot  solidi¬ 
fication  was  determined  by  coarsening  ("ripening").  A  simple  analysis, 
based  on  coarsening  kinetics,  was  found  to  be  of  engineering  applic¬ 
ability  to  predict  effects  of  alloy  analysis  and  other  solidification 
variables  on  dendrite  arm  spacing^. 

During  the  period  for  which  this  report  is  written,  research  has 
been  continued  on  effects  of  coarsening  kinetics  on  dendrite  structure, 
and  on  structure  in  rapidly  solidified  alloys.  In  addition,  limited 
work  has  been  conducted  on  solidification  in  composites,  and  on  effects 
of  melt  agitation  on  structure.  Finally,  a  significant  effort  has  been 
initiated  to  study  mechanical  properties  of  ultra-high  strength  wrought 
aluminum-base  alloys,  produced  using  solidification  techniques  to 
assure  complete  homogeneity  of  the  original  ingot.  These  studies  are 
to  be  summarized  in  a  later  report;  this  report  will  deal  with  work 
that  has  comprised  the  major  effort  during  this  contract  year  -  i.e., 
macrosegregation  in  casting  and  ingot  solidification. 
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This  work  has  provided,  for  the  first  time,  a  basis  for 
understanding  a  wide  variety  of  apparently  different  types  of 
macrosegregation.  The  basic  analysis  presented  provides  quantitative 
as  well  as  qualitative  understanding  of  effects  such  as  inverse 
segregation,  centerline  segregation,  and  banding,  The  results  should 
be  of  considerable  engineering  usefulness  in  understanding  and 
controlling  macrosegregation  in  ferrous,  as  well  as  non-ferrous  alloys. 

Theoretical  portions  of  this  work  employ  a  model  for  solute 

redistribution  in  dendritic  solidification  discussed  in  detail  by 

13  2  4 

Brody  and  Flemings  *  ,  and  Bower,  Brody,  and  Flemings  *  .  Concepts 

employed  in  considering  the  interdendritic  fluid  flow  are  those 

described  by  Piwonka  and  Flemings^.  Theoretical  work  of  Schell** 

12 

and  Kirkaldy  and  Youdelis  on  inverse  segregation  has  served  as  an 
important  foundation  for  thus  study.  Finally,  special  acknowledgement 
is  made  of  stimulating  discussions  with  Professor  John  Cahn  of 
Massachusetts  Institute  of  Technology. 
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I.  INTRODUCTION 


In  the  wtrks  of  Vannoccio  Biringuccio^  published  in  1540,  one 
of  the  first  references  to  macrosegregation  is  noted.  He  refers  to 
the  problem  of  exudation  and  inverse  segregation  encountered  in  the 
manufacture  of  bronze  gun  barrels.  From  that  day  to  this,  the  pro¬ 
blem  of  macrosegregation  has  been  studied  by  many  researchers.  The 
usual  approach  has  been  to  isolate  one  type  of  macrosegregation 
(l.e.  centerline  segregation,  negative  cone  of  segregation,  exudation, 
inverse  segregation,  banding,  etc.)  with  the  exclusion  of  all  other 
types.  The  result  has  been  a  multiplicity  of  theories,  each  one 
designed  to  explain  a  specific  type  of  macrosegregation.  This  section 
will  consider  many  of  these  macrosegregation  effects,  present  some  of 
the  pertinent  literature  and  discuss  theories  pertaining  to  each 
effect.  In  a  subsequent  section  it  will  be  shown  that  many  of  the 
different  types  of  macrosegregation  are  not  due  to  unrelated  phenomena 
but  can  be  understood  simply  in  terms  of  the  direction  and  velocity  of 
solute  rich  liquid  which  flows  to  feed  solidification  and  thermal 
contractions. 

A.  Inverse  Segregation 

One  of  the  most  thoroughly  investigated  macrosegregation  effects 
is  inverse  segregation.  In  the  late  1920’s  and  early  1930's  many 

theories  were  proposed  to  explain  the  mechanism  of  inverse  segregation. 
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Brick  and  Phillip  have  reviewed  these  proposals  and  rejected  most  of 
the  suggested  mechanisms  as  being  either  speculative  or  unconfirmed  by 
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experimental  observations.  They  did  accept,  as  most  observers  do 
today,  that  Inverse  segregation  Is  due  to  solidification  and  thermal 
volume  contraction  with  resulting  flowback  of  solute  enriched 
liquid,  increasing  local  average  composition  above  Cq  (average  melt 
composition) . 

Most  of  the  definitive  work  which  substantiates  the  volume 

contraction  and  interdendritic  flow  theory  for  explaining  inverse 
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segregation  was  published  in  the  1940' s.  Adams  studied  inverse 

segregation  in  aluminum-copper  alloys  (3  per  cent  and  7  per  cent 

copper)  and  reported  extensive  experimental  results.  However,  he, 
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along  with  Sauerwald  and  Masing  ,  while  attempting  quantitative 
predictions,  derived  no  analytical  expressions.  The  first  analytical 
work  to  predict  inverse  segregation  which  could  be  checked  experiment¬ 
ally  was  proposed  by  Schell^.  His  expression  predicted  maximum 
segregation  at  the  chill  face  as  a  function  of  alloy  content.  Good 
order  agreement  was  obtained  with  experiment.  Kirkaldy  and  Youdelis^ 
recognized  that  the  basic  mechanism  that  explains  inverse  segregation 
accounts  for  the  details  of  the  segregation  distribution  in  a  uni- 
directionally  solidified  ingot.  They  correlated  previously  available 
data  on  inverse  segregation  in  aluminum-copper  alloys,  and  developed 
expressions  in  differential  form  to  predict  macrosegregation  throughout 
the  length  of  a  unidirectionally  solidified  ingot;  these  were  solved 

g 

by  numerical  analysis.  Later  Youdelis  and  Colton  extended  the  work  to 
a  study  of  inverse  segregation  in  aluminum-zinc  Ingots.  A  limiting 
case  of  the  macrosegregation  expression,  given  later,  is  exactly  that 


3 


of  Klrkaldy  and  Youdelis^;  in  addition,  it  is  shown  that  by  making 
certain  approximations,  a  simple  analytical  solution  is  possible. 

The  literature  abounds  with  investigations  on  inverse 
segregation.  Some  of  the  data  are  of  very  early  origin  and  experiment¬ 
al  procedures  and  results  are  questionable.  Extensive  reviews  of  the 
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literature  ate  also  available,  notably  by  Vaughn  reviewing  work  to 
1937,  by  Pell-Walpole^  covering  work  to  1949  and  most  recently  by 
Vosskuhler^*  where  65  references  are  cited,  many  recent,  on  inverse 
segregation.  In  the  theoretical  section  it  will  be  shown  how  not 
only  inverse  segregation  but  also  the  "negative  cone"  of  segregation, 
banding,  centerline  segregation,  and  localized  high  or  low  solute 
regions  can  be  related  to  the  same  mechanism  which  causes  inverse 
segregation;  namely  flow  of  solute  enriched  liquid  to  feed  solidifica¬ 
tion  and  thermal  contraction.  For  this  section,  however,  existing 
theories  on  these  effects  will  be  considered. 

B.  Negative  Cone  of  Segregation 

Probably  no  macrosegregation  phenomena  is  less  understood  than 

the  cone  of  negative  segregation  which  often  appears  at  the  bottom 

of  large  ingots.  The  usual  explanation  for  this  macrosegregation 

effect  is  related  to  initial  crystals  which  form  being  of  lower 

composition  than  the  average  composition  of  the  melt.  As  the  crystals 

form  they  are  said  to  settle  to  the  bottom  of  the  ingot  (since  they 

are  more  dense  than  the  liquid)  and  pile  up  in  a  cone  of  high  purity, 

12 

low  solute  composition  .  Several  difficulties  exist  with  the  theory: 


1.  In  an  ingot,  the  llquidus  temperature  is  reached  by  most  of 
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the  bulk  liquid  before  solidification  starts  * 

2.  It  can  be  shown  from  Stokes'  law  that  the  rate  of  descent  of 
crystals,  in  liquid  metals,  compared  to  solidification  times, 

is  too  slow  to  descend  the  required  distance  in  the  available 

.  *  15 

time  . 

One  of  the  new  ideas  on  this  phenomena  is  reported  by  Marburg^. 
He  associates  the  "negative  cone"  with  accelerated  vertical  and  trans 
verse  solidification  after  the  initiation  of  solidification.  In  a 
subsequent  section  this  will  be  discussed  and  related  to  the  flow  of 
solute  rich  liquid  to  feed  solidification  contraction. 

C.  Centerline  Segregation 

Another  macrosegregation  effect  little  understood  is  centerline 
segregation.  While  data  are  available  on  the  magnitude  of  the  segre¬ 
gation  encountered^,  theories  usually  are  limited  or  speculative. 
Until  recently,  the  generally  accepted  theory  related  centerline 

segregation  to  a  buildup  of  solute  ahead  of  the  interface  of  the 
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columnar  grains  growing  away  from  the  vertical  wall  '  .  As  the 

grains  grow  toward  the  center,  the  rejected  solute  was  visualized  as 

being  pushed  ahead  resulting  in  a  solute  rich  centerline  where  the 

growth  fronts  converge.  Recently,  however,  it  has  been  shown  that 

there  can  be  no  significant  solute  buildup  in  front  of  growing  tips 

of  dendrites  at  the  growth  rates  and  gradients  encountered  in  usual 
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casting  processes  '  *  .  In  a  following  section  the  mechanism  for 


5 


centerline  segregation  will  also  be  considered  a?  being  related  to 
solute  rich  fluid  flow  to  feed  solidification  contraction. 

D.  Banding 

Probably  no  problem  gives  the  ingot  maker  or  foundryman  more 

problem  than  banding  type  macrosegregation.  Banding  refers  to  regions 

where  an  Impurity  of  solute  in  the  base  alloy  is  localized  in  a  narrow 

band  or  strip.  In  ferrous  alloys  these  bands  may  be  high  in  Inclusion 

content,  but  in  all  alloys  are  of  somewhat  different  composition  than 

the  base  alloy  (typically  .05  to  1.00  per  cent,  depending  on  base 

solute).  Banding  occurs  in  large  ingots,  centrifugal  castings  and 

continuously  cast  billets,  and  often  seems  to  occur  in  a  random 

manner.  In  large  steel  ingots,  Marburg^  has  associated  inverted  V 

and  V-type  banding  segregation  with  changes  in  the  rate  of  solidifica- 
21 

tion.  Takenaka  shows  that  banding  type  defects  in  centrifugally 

cast  steel  pipe  are  influenced  by  vibration  in  the  centrifugal 

casting  system  and  solidification  rate  in  the  melt.  This  appears  to 

be  as  close  as  anyone  comes  to  relating  banding  to  a  mechanism  similar 

22 

to  that  of  inverse  segregation.  Northcott  suggests  that  vibration 

might  cause  banding  by  affecting  the  boundary  layer  in  front  of  the 

growing  dendrites.  In  view  of  the  very  small  thickness  of  this 
20  37  38 

layer  *  *  ,  however,  it  seems  unlikely  that  this  theory  is  correct. 

23 

Cumberland  proposes  a  theory  based  on  the  bands  being  incipient 
laps  resulting  from  irregularities  in  the  manner  of  liquid  metal 
flow  within  the  rotating  mold.  This  theory,  however,  again  requires 
extensive  buildup  of  solute  at  the  tips  of  the  solidifying  interface. 
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An  older  theory  by  Howson4,  relates  banding  to  a  mechanical 
separation  of  elements  due  to  centrifugal  forces.  The  theory, 
however,  does  not  take  Into  account  solidification  time  which  can 
be  quite  rapid. 

Finally,  while  it  is  not  strictly  a  banding  effect,  consideration 

will  be  given  to  the  segregation  in  continuously  cast  billets  and 
25 

slabs.  Wolf  discusses  this  problem  in  Al-Zn-Mg-Cu  alloys  and 
reports  on  similar  investigations  by  other  researchers.  While  much 
data  are  reported  and  conclusions  drawn  as  to  what  affects  the 
segregation  in  continuous  castings,  little  seems  to  be  available 
concerning  possible  mechanisms.  Compositions  toward  the  center  of 
continuously  cast  ingots  were  found  to  decrease  which  can  again  be 
related  to  solute  flow  to  feed  shrinkage. 

E.  Solute  Enriched  Regions 

A  type  of  macrosegregation  related  to  banding  is  the  phenomenon 
of  localized  high  solute  regions  usually  found  in  castings  and  ingots. 
By  localized  enrichment  is  meant  any  region  where  solute  content  can 
be  significantly  higher  than  average  base  alloy  composition.  The 
types  of  macrosegregation  effects  usually  associated  with  this 
phenomena  are  exudation,  hot  tearing  type  defects  and  solute  enrich¬ 
ment  under  risers.  These  three  phenomena  will  be  considered 
separately  even  though  they  are  associated  effects,  related  by  the 


mechanism  involved. 


7 


1.  Exudation. 

Exudation  is  usually  found  in  conjunction  with  Inverse 

segregation,  and  the  two  effects  are  often  confused.  Exudation 

occurs  at  the  mold-metal  or  chill-metal  interface,  apparently  as  a 

result  of  metal-mold  gap  formation.  After  solidification  is  started 

at  a  wall  or  chill,  contraction  of  the  metal  causes  an  air  gap  to 

form.  When  solute  enriched  liquid  is  forced  from  between  dendrites 

to  fill  the  gap,  a  layer  of  enriched  solid  which  can  approach 

eutectic  composition  is  formed.  Exudation  has  been  studied  and  th< 

26  27 

variables  discussed  which  affect  thickness  of  the  layer  *  and  its 
28  29 

formation  *  .  All  data  indicate,  that  exudation  is  related,  as 

are  many  macrosegregation  effects,  to  the  flow  of  solute  enriched 
interdendritic  liquid. 

2.  Solute  Enrichment  Under  Risers. 

One  of  the  effects  readily  related  to  solute  enriched  fluid  flow 

is  macrosegregation  beneath  casting  risers  and  ingot  "hot  tops". 

At  the  later  stages  of  solidification  (i.e.,  when  casting  or  ingot  is 

almost  completely  solid)  some  solidification  has  already  occurred 

in  the  riser  or  hot  top.  Because  of  the  partial  solidification  in 

the  riser  or  hot  top  the  liquid  has  become  enriched  in  solute.  When 

further  contraction  occurs  in  the  casting  or  ingot  the  solute 

enriched  liquid  feeds  from  the  riser  to  the  casting.  The  result  is 

the  greatest  enrichment  in  the  casting  or  ingot  at  the  last  location 

30 

to  solidify  (directly  under  the  hot  top  or  riser).  Bishop  et  al 
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have  reported  extensively  on  this  effect  In  steel  castings.  It 

Is  presently  being  studied  at  the  MIT  casting  and  solidification 
31 

section  using  aluminum-4. 5  per  cent  copper  alloy.  It  will  be 
seen  later  in  this  report  that  this  under-riser  segregation  is 
expected  to  be  most  severe  when  the  riser  is  of  inadequate  size, 
or  inadequately  insulated. 

3.  Hot  Tearing. 

The  final  type  of  macrosegregation  to  be  considered  is  "filled 
hot  tears".  "Hot  tearing"  is  associated  with  the  solid  contraction 
which  occurs  when  a  casting  or  ingot  cools.  When  the  contraction  is 
hindered  by  the  mold,  cores,  or  differences  in  casting  section  size, 
complete  rupture  of  a  segment  of  the  casting  or  ingot  may  occur. 

When  these  fractures  occur  within  the  liquid-solid  zone,  they  are 
termed  "hot  tears"  and  when  an  incipient  tear  is  fed  from  a  surround¬ 
ing  region,  the  result  is  a  "filled  hot  tear",  generally  rich  in 

32 

solute.  Rosenberg  et  al  have  studied  hot  tearing  in  aluminum- 

copper  alloys  and  report  data  which  can  be  accounted  for  by  solute 

rich  flow  feeding  the  "hot  tear".  While  much  work  has  been  done  on 
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the  fundamental  mechanism  of  "hot  tearing"  in  steelJ  '  and  non- 
35  36 

ferrous  alloys  *  little  attention  has  been  given  to  the  resulting 
solute  enriched  flow  necessary  when  the  "hot  tear"  is  fed. 


In  summary,  five  different  types  of  macrosegregation  with  varia¬ 
tions  have  been  discussed,  primarily  on  the  basis  of  generally  accepted 
exsiting  theories.  In  the  next  section  it  will  be  shown  that  all  these 
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effects  are  related  to  solute  enriched  fluid  flow  feeding 
solidification  and  thermal  contraction,  and  that  they  can  be 
quantitatively,  as  well  as  qualitatively,  understood  on  this  basis. 


10 


II.  THEORETICAL  DEVELOPMENT 

37  38 

In  two  recent  papers  '  ,  the  problem  of  microsegregation  in 

dendritic  solidification  was  considered.  It  was  concluded  that  for 
usual  casting  and  Ingot  solidification,  microsegregation  in  binary 
alloys  can  accurately  be  predicted  by  the  "non-equilibrium  freezing 
equation"  provided  there  is  no  diffusion  in  the  solid  and  no  macro¬ 
segregation.  Up  to  the  limit  of  solubility  of  the  primary  solid 
phase  this  is  written: 


(l-k)CLdfL  -  -fLdcL 


(1) 


where:  *  liquid  composition  within  the  "volume  element" 

(a  small  volume  the  order  of  size  of  the 
dendrite  arm  spacing) 


“  weight  fraction  liquid  in  the  volume  element 
k  ■  equilibrium  partition  ratio 


For  the  special  case  of  constant  partition  ratio,  equation  (1) 


is : 


(2) 


or,  since  C*  ■  kC^  and  fg  ■  (1  -  f^) 


C*  -  kC0(l  -  fs) 


(3) 
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where:  C*  ■  solid  composition  at  the  liquid-solid  Interface 
fg  ■  weight  fraction  solid 

Equation  (3),  for  the  case  of  no  solid  diffusion,  describes  solid 

composition  as  a  function  of  weight  fraction  solid  after,  as  well  as 

during,  solidification.  It  therefore  completely  describes  "severity 
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of  microsegregation",  although  not  its  geometry  * 

Here  consideration  is  first  given  to  how  the  foregoing,  and 
particularly  equation  (1),  is  modified  by  mass  flow  in  or  out  of  the 
volume  element  during  solidification.  The  net  change  in  composition 
of  the  element  at  the  end  of  solidification  is  then  calculated  (l.e., 
"macrosegregation") . 

PART  I:  GENERAL  MACROSEGREGATION  EXPRESSIONS 

A.  Local  Solute  Redistribution 

Consider  a  small  volume  element  at  coordinates,  x»  y,  z  within  a 
solidifying  ingot.  These  coordinates  are  to  be  such  that  no  solid 
material  enters  or  leaves  the  volume  element  during  solidification. 
Hence  tie  location  x,  y,  z  may  move  with  respect  to  the  mold  wall  as 
a  result  of  thermal  contraction  or  other  solid  movement,  but  refers 
to  the  same  local  region  of  a  given  dendrite.  Solute  enters  or 
leaves  the  element  only  by  liquid  flow  to  feed  shrinkage.  Mass  flow 
by  diffusion  is  neglected;  such  mass  flow  has  been  shown  to  be 
negligible  in  the  absence  of  convection  sufficiently  vigorous  to 
penetrate  interdendritlc  regions.  Local  perturbations  of  the  macro¬ 
scopic  composition  gradient  are  assumed  negligible  so  C.  is  uniform 
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within  the  volume  element.  This  model  Is  Identical  to  that  described 
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earlier  ’  ,  except  that  solidification  contraction  and  thermal 

contraction  are  considered.  The  change  In  total  solute  In  the  element 
per  unit  time  during  solidification  Is  then: 

dT  (p£)  *  ",'bl8lci?  <4) 

(a)  (b) 

3 

where:  p  *  local  average  density  (of  liquid  and  solid),  gm/cm 

C  ■  local  average  composition  of  liquid  and  solid, 
gm/gm 

3 

*  liquid  density,  gm/cm 

v  *  local  velocity  of  liquid  relative  to  solid  (resulting 
from  solidification  contraction) ,  cm/sec 

g^  ■  volume  fraction  liquid 
The  two  terms  are: 

(a)  Total  local  rate  of  change  of  solute  during  solidification 

(b)  Local  rate  of  change  of  solute  resulting  from  liquid  flow 
relative  to  the  solid. 

The  change  in  total  mass  in  the  element  per  unit  time  is: 


dt 


-7-P[gLv 


(5) 


and  this  change  is  the  sum  of  volume  changes  in  the  liquid  and  solid 


during  solidification: 
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dt  “  dt  ^PS8S  +  pLgL^ 


(6) 


where:  Pg  ■  solid  density ,  gms/cm' 
gg  ■  volume  fraction  solid 


In  the  absence  of  pore  formation: 


8S  +  gL  *  1 


(7) 


The  term  on  the  right  of  equation  (4)  may  be  rewritten: 


"V  *  PL8LCLV  "  “  CLV  *  pLRLV  "  pl8LvVCL 


(8) 


Substituting  equations  (3)  -  (8)  in  (4): 


dt  ^SPS8S  +  CLPL8!^ 


CL  dt  ^PS8S  +  PL8L*  ”  PL8LV7CL  ^ 


Assuming  the  partition  ratio  applies  at  the  Interface,  no 
diffusion  in  the  solid,  and  constant  solid  density,  Pg: 


d^SpS8S^  “  kCLpSd8S 


(10) 


Combining  (9)  and  (10): 


dgI 

pscl(1  -  k)  dT 


dCL 

~PLgL  ^  dt  * 


pl8lvVCl 


(11) 
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or: 


VfL 

gLdt 


where  6  is  the  "solidification  shrinkage": 


(12) 


Equation  (12)  is  the  basic  local  solute  redistribution  equation 
used  hereafter  to  calculate  macrosegregation.  It  is  written  for 
the  general  case  of  three  dimensional  heat  and  fluid  flow,  for  (a) 
constant  solid  (but  not  liquid)  density  during  solidification,  and 
(b)  negligible  net  solute  change  from  diffusion. 


B.  Local  Solute  Redistribution.  Planar  Isotherms 


Assume  isotherms  are  planar  and  perpendicular  to  x  axis  where 
x  is  the  distance  from  the  mold  wall.  Liquid  composition  at  the 
liquid-solid  interface  is  a  function  only  of  temperature  provided 
undercooling  from  solute,  radius  of  curvature,  and  kinetics  effects 
is  negligible.  Hence,  planar  isotherms  require  planar  isoconcen¬ 
trates  and: 


i  ( 


3C, 

dX 


(13) 


also: 


-  iv  +  jv  +  kv 
x  y  z 


v 


(14) 
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Combining  (13)  and  (14)  with  (12): 


dC. 


(  1  -  k  X  CLdgL  . 

“  V  1  _  Q  /  ^  “  v 


1  -  6 


St 


3C, 

9x 


)  v  dt 
x 


(15) 


Now 


is  a  function  only  of  temperature.  Isotherms  and 


dx 


laoconcentrates  move  with  velocity  “  at  x.  Hence, 


and: 


3C_  j 
/  L  v  dfoc 

K  9>:  '  dt 


fi 

dt 


(16) 


9x 


dC, 

dx 


(17) 


Similarly: 


il 

dx 


dT 

dx 


(18) 


where:  T 


temperature 


9T  dT  dT 

Now,  let  G  ■  ”  dx  (temPerature  gradient),  and  e  ■  —  (rate  of 

temperature  change),  then: 


_e 

G 


dx 

dt 


(19) 


Combining  (17),  (18),  and  (15): 


-  ( 


1  -  k 


v  G  dC. 

><1  +  -r>  r 


(20) 


Equation  (20)  is  the  solute  redistribution  expression  (equation  12) 
for  the  case  of  unidirectional  heat  flow,  hut  where  fluid  flow  is  not 
necessarily  unidirectional,  k  and  6  are,  in  general,  functions  of  C.  . 
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G  and  c  (at  x,  y,  z)  are  also  functions  only  of  temperature,  and 

hence  of  C^;  these  quantities  are  readily  measured  experimentally,  or 

In  principle,  calculated  from  heat  flow  considerations.  Equation  (20) 

becomes  identical  to  that  of  Klrkaldy  and  Youdelis^  when  v  ■  v  "0. 

y  z 

Evaluation  of  v  is  discussed  below, 
x 

C.  Evaluation  of  Flow  Velocity  Vector.  Planar  Isotherms 
1.  Flow  Perpendicular  to  Isotherms. 


Consider  first  the  case  of  flow  only  perpendicular  to  planar 
isotherms.  Equation  (5)  and  (6)  then  become,  for  locations  x' : 


£ 

dt 


(pXv;> 


Tt  <ps*s  +  PW 


(21) 


where:  x'  is  distance  from  the  chill  face  to  the  location 
considered. 


When  the  chill  face  is  above  the  non-equilibrium  solidus,  T  , 

0  <  x'  <  x,  and  when  the  chill  face  is  below  T^.,  x^,  <  x'  <  x.  The 

primes  on  other  quantities  denote  that  these  are  at  x'.  Now 

assuming  constant  p,,  and  considering  the  case  where  the  chill  face 

(x'  ■  0)  is  above  Tr: 

t 


°LgLVx 


/  3(PL8Lvx) 
0 


x'-x 

1  K 

x'-O 


dgL  X  "X  dpL 

°S)  dt  3x'  +  !  gL(  dT  )  9x' 

x'-O 


or 


PL8LVx 


X  -X 

/ 

x'-O 


K 

<ps-°[)(dr)(3x,) 


x'“x  dp' 

/  g‘(  ax' 

x'-O 


(22a) 
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when  the  chill  face  is  below  the  non-equlllbrlum  solidus,  In  the 
case  of  an  alloy  containing  some  eutectic,  eutectic  of  volume 
fraction  gp  precipitates  at  the  non-equilibrium  solidus,  Tp.  Then: 


pLgLVx 


'  (ps-pi)(  df)  3x’ 

X'-Xp 


dpL 

(  — -  ) 
k  dt  1 


3x'  - 


^PSE  "  PLE^  8EUE  (22b) 

where:  p,,  ,  Ppp  ■  densities  of  solid  and  liquid  eutectic 

respectively 

U  *  velocity  of  eutectic  isotherm,  T 
E  E 

dg.'  is  a  function  of  x'  as  well  as  macrosegregation  at  x'. 

L 

Hence,  to  rigorously  determine  local  solute  redistribution  at  x, 
equation  (22)  and  equation  (30)  must  be  solved  by  a  tedious  iterative 
procedure.  To  permit  simpler,  analytical,  solutions  for  the  segrega¬ 
tion,  an  alternate  procedure  is  employed  in  much  of  the  following. 

Pg,  Pp  and  k  are  assumed  constant,  and  for  velocity  calculation,  g^ 
is  assumed  to  vary  linearly  with  temperature  within  the  liquid-solid 
zone. 


1.  Flow  Not  Perpendicular  to  Isotherms. 

Calculation  of  the  general  case  of  macrosegregation  when  fluid 
flow  is  not  perpendicular  to  isotherms  requires  detailed  knowledge 
of  v^  ■  f(x',  y',  ::’).  In  the  special  case  where  fluid  flow  is  not 
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unidirectional  but  where  v'  and  g'  are  independent  of  y'  and  z',  an 

X  L 

expression  similar  to  equation  (22b)  is  written  (for  the  chill  face 

below  T_) : 

E 


APLgLV* 


/  (P' 


P[>(  df> 


■i  <  a r  )A' *’ 


iC  -  Cs  -  C0  (24) 

where:  Cg  ■  average  composition  at  x,  y,  z  after  solidification 
Cq  *  bulk  composition 

"Positive  segregation"  refers  to  Cg  >  CQ,  and  "negative 
segregation"  for  Cg  <  C^.  For  a  binary  alloy  containing  a  primary 
phase  a  and  some  eutectic,  the  segregation  after  solidification  is 
given  by: 


where:  p^,  pg^  are  densities  of  a  and  eutectic,  respectively 

(p  ■  p_  ■  constant) 
a  b 

C  ■  average  composition  of  the  entire  amount 

of  primary  a  phase  in  the  final  solid  ingot 

Cg  ■  composition  of  eutectic 


The  average  composition,  Cg,  of  any  portion  of  the  a  phase,  gg, 
is  given  by. 

8S  +  A8S 

PS^SAgS  “  PS  ^  CSdgS  <26) 

8S 

and  the  mass  of  solute  in  the  entire  primary  phase,  Pg8aCa  determined 
by  summing  equation  (26)  for  all  portions  of  the  phase  Agg,  which 
total  to  g  .  Cc  as  a  function  of  gc  is  determined  from  equation  (20) 

Ol  b  b 

recalling  that  Cg  ■  kCL,  and  gg  ■  1  -  g. . 


Incomplete  Feeding 


In  the  foregoing  and  in  most  of  the  following,  it  is  assumed  that 
feeding  is  complete.  However,  under  many,  if  not  most,  solidification 
conditions,  pores  are  expected  to  develop  in  the  late  stages  of 
solidification  as  a  result  of  resistance  to  fluid  flow,  precipitation 

OQ 

of  dissolved  gas  or  a  combination  of  both  .  The  preceedlng  analyses 
are  readily  adaptable  to  treatment  of  macrosegregation  under 


these  conditions. 


20 


In  the  simplest  case,  feeding  Is  assumed  to  stop  completely  at 
a  critical  fraction,  liquid,  g^»  at  time  tv,  location  x^  and  liquid 
composition  becomes  C^.  Hence  flow  velocity  to  feed  shrinkage  will 
be  less  In  the  case  of  pore  formation,  since  the  lower  Integration 
limit  on,  for  example,  equation  (22b)  becomes  x  Instead  of  x  . 

Further,  In  evaluation  of  the  solute  redistribution  expression, 
equation  (20),  one  set  of  limits  now  becomes  g^  and  C^;  for 

y 

g^  <  no  further  feeding  occurs.  Final  macrosegregation  is  then 
simply  calculated  from  summing  the  solute  In  the  solid  and  liquid 
phases  at  tv  (since  after  t^  no  further  solute  enters  or  leaves  the 
volume  element,  although  solidification  progresses). 

PART  II:  UNIDIRECTIONAL  SOLIDIFICATION 

As  examples  of  application  of  the  foregoing  to  calculations  of 
macrosegregation,  the  case  of  "unidirectional"  solidification  is 
first  considered;  i.e.,  solidification  with  planar  isotherms, 
unidirectional  heat  flow.  Some  limiting  cases  are  first  considered, 
and  then  segregation  versus  distance  from  the  chill  wall  Is  examined 
for  various  conditions  of  heat  and  fluid  flow. 

A.  Limiting  Cases 

As  examples  of  the  foregoing,  solutions  to  equation  (2Q)  for  three 

limiting  cases  are  considered  first  here:  (a)  vxG/e  ■  0,  (b)  "steady 

state  solidification",  and  (c)  non  "steady  state",  but  finite  vxG/e 

4  f(gT).  8  and  8  are  assumed  constant  and  equal. 

L  6 
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1. 


0. 


At  a  chill  face,  v  ■  0,  and  provided  G/e  is  finite,  v  G/e  ■  0. 

X  X 

For  constant,  k,  6  and  no  solid  movement  during  solidification,  equa 
tlon  (20)  becomes: 


and 


/ 

1 


1  -  B 
1  -  k 


(27) 


(28) 


Here,  increasing  B  results  in  increasing  at  a  given  g^  and 
hence  in  greater  positive  macrosegregation. 

2.  Steady  State  Solidification. 

Consider  now  unidirectional  "steady  state"  solidification  of  a 
long  ingot;  i.e.,  planar  isotherms  move  at  a  constant  velocity 
U  *  -  c/G,  along  the  x'  direction  in  an  ingot  of  constant  cross 
section,  and  chill  face  is  below  the  non-equilibrium  solidus.  For 
B  ■  B  *  constant,  equation  (22b)  reduces  to: 


Substitution  of  equations  (29)  in  (20)  and  integrating  from 
the  limit  gL  ■  1  at  ■  CQ  to  g^,  yields: 


22 


g 


L 


(30) 


Equation  (30),  for  the  special  case  of  0  ■  B  ■  constant  yields, 
at  the  end  of  solidification,  precisely  the  same  solute  distribution 
as  equation  (2).  The  final  composition  (weight  fraction  or  volume 
fraction  solute)  is  exactly  that  of  the  original  liquid  and  no  macro¬ 
segregation  results. 

v  G  v  G 

3.  Finite  —  ,  Non  "Steady  State",  4  f(gT). 

C  C  L 

v  G 

x 

As  will  be  shown  later  — —  can  vary  substantially  from  that 

value  at  the  chill  face  (0  for  finite  G/e)  ,  or  that  value  at  steady 

6  e  VXG 

state  (  - t  )(  —  ).  Positive  values  of  —  ,  different  than  the 

1  •  P  b  E 

steady  state  value  result  from  (a)  cross  sectional  area  change  of 

the  ingot,  (b)  change  in  solidification  rate  at  locations  between  the 

non-equilibrium  solidus  and  x,  (c)  change  in  G  or  c  (at  x) .  Negative 
v  G 

X 

values  of  •  ~  ~~  result  when  liquid  feed  metal  flows  in  the  same  direc¬ 
tion  as  the  direction  of  movement  of  isotherms  at  x,  as  by  reheating 
(remelting)  at  locations  between  the  non-equilibrium  solidus  and  x. 
Examples  other  than  those  mentioned  above  will  be  considered  later. 


v  G 

For  the  special  cases  where  -  4  f(gT)»  and  where  the  chill 

E  L 

face  is  below  the  solidus  temperature,  equation  (20)  can  be  directly 


integrated: 


g 


L 


1  -  B 
1  -  k 


v  G 

)(  1  +  “ ~  ) 


(31) 
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U .  Limiting  Cases.  Numerical  Examples 

Consider  for  numerical  example  of  the  foregoing  aluminum-4.5 
per  cent  copper  alloy,  solidifying  unidirectlonally  from  a  chill  face. 
The  appropriate  portion  of  the  phase  diagram  is  shown  in  Figure  la; 
Figure  lb  shows  schematically  the  liquid-solid  "mushy"  zone  during 
solidification.  Dendrite  tips  are  at  the  liquldus  temperature,  T^, 
and  liquid  exists  within  the  dendrite  arms  to  the  eutectic  tempera¬ 
ture  Tg  (Figure  lc).  It  should  be  noted  that  while  dendrites  are 
shown  schematically,  the  following  analysis  is  independent  of 
dendrite  morphology. 

From  the  phase  diagram,  k  *  0.172,  maximum  solid  solubility  at 
the  eutectic  is  5.65  per  cent  copper,  and  eutectic  composition  is  33 
per  cent  copper,  and  assume  6  *  $e  ■  .055,  a  constant.  (6  is  not,  in 
fact,  a  constant  for  this  alloy;  however,  the  assumed  value  of  .055 
results  in  calculated  segregation  in  good  agreement  with  that  observed 
experimentally,  as  will  be  shown  later.) 

Using  the  foregoing  data,  the  solute  redistribution  equations 

for  each  of  the  foregoing  limiting  cases  are  plotted  in  Figure  2; 

10 B 

that  for  the  third  case  is  plotted  for -  ■  - —  .  In  each  case  the 

G  1  —  P 

final  average  composition  at  the  location  in  question  is  given  by  the 

area  under  the  curve  drawn;  this  average  composition  decreases  with 
v  G 

increasing  — - —  . 

This  final  average  composition  is  readily  determined  graphically 
or  by  solution  of  equation  (26).  Figure  3  shows  the  general  result 


ATOMIC  PERCENT  COPPER 
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Figure  la:  Aluminum  rich  portion  of  aluminun-copoer  oh 


Figure  lb:  Schematic  diagram  of  dendrites  in  liquid-solid  mushy  zone 
during  solidification. 


Figure  lc:  Temperature  distribution  in  unidirectional  solidification. 


COMPOSITION  SOLID  (PERCENT  COPPER) 
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Figure  2:  Solute  redistribution  for  three  solidification  cases 
(a)  vxG/e  ■  0,  i.e.,  chill  face  solidification, 

(h)  vxG/e  *  8/(1  -  8),  i.e.,  steady  state  solidification, 
(c)  vxG/t  ■  108/(1  -  8),  i.e.,  10/1  change  in  solidifica¬ 
tion  cross  section. 
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of  average  composition  C_  (and  segregation,  AC)  plotted  against  the 

VXG  i  _  c 

dimensionless  group  of  variables,  p  ■  - -  ■  ■■  .  For  the  case  of 

e  P 

the  chill  face,  p  ■  0,  for  steady  state,  p  ■  1.  In  the  general 
case,  however,  p  can  have  a  wide  range  of  positive  or  negative 
values  with  the  marked  segregation  Illustrated. 

C.  Macrosegregation  Versus  Distance  from  Chill 
1.  Description  of  Model,  Constant  k,  8. 

We  consider  first  the  case  of  macrosegregation  versus  distance 
from  the  chill,  assuming  constant  k  and  8.  Assume,  as  in  Figure  4a, 
liquidus  and  non-equilibrium  solidus  (eutectic)  isotherms  move 
according  to  the  general  relations: 

\  ‘  +  aL  <32a) 

Xg  «  nEtb  +  *E  (32b) 

where  x.  and  x_  are  the  distance  from  the  chill  to  the  liquidus 
and  eutectic  isotherm  respectively,  t  is  time  and  constants  are  as 
defined  by  Figure  4a.  By  proper  choice  of  a^  and  a^.  we  can  readily 
use  relations  derived  above  for  a  variety  of  unidirectional  solidi¬ 
fication  cases  including  (1)  solidification,  at  constant  rate 

1/2 

and  gradient  (Figure  4b),  (2)  xT  and  x.  proportional  to  t  as 

L  L 

expected  for  unidirectional  solidification  against  a  chill  wall  with 

no  mold-metal  interface  resistance  (4c)  and  (3)  x.  and  x_  propor- 

L  L 

1/2 

tlonal  to  t  plus  a  constant;  expected  for  unidirectional  solidifi- 

38 

cation  against  a  chill  wall  with  mold-metal  interface  resistance 


t  t 

/  ✓ 

(TIME  )b 

'  °L 
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L  ■  f 
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/ 

* 

(a) 

Figure  4(a) : 

Liquidus  and 

eutectic  isotherm  movement ,  interface 

resistance  case,  n  4  nr. 

Li  L 
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T1IHD  IMOdd  30NV1SI0 


“IUHO  IMOdd  3DNV1SIQ 


"HIHD  IMOdd  3DNV1SIQ 


Figure  4:  Liquidus  and  eutectic  isotherm  movement:  (b)  steady  state,  nE  =  n^.  =  constant 
(c)  no  interface  resistance,  a£  =  a£  =  0,  (d)  general  case,  aL  4  a|,  nL  4  nE. 
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For  analytic  solution,  it  is  now  necessary  to  assume  a 
distribution  of  fraction  liquid,  versus  x'  in  the  mushy  zone. 
Reasonable  alternatives  include: 


(1)  A  distribution  as  given  by  the  steady  state  solute 
redistribution  expression,  equation  (30).  Such  a  distribution 
would  clearly  be  correct  for  linear  temperature  distribution 
and  small  macrosegregation.  This  is  sketched  in  Figure  5a. 

(2)  Linear  gl  versus  x',  assuming  g  constant  and  given  by 

L  t 

equation  (30) (g  •  .0902  for  aluminum-4.5  per  cent  copper), 
Figure  5b. 

(3)  Linear  gT'  versus  x’,  assuming  g  *  0  at  x_.  Figure  5c. 

L  L  L 


An  analytic  solution  for  segregation  based  on  condition  (1)  above 

would  be  quite  complicated.  A  solution  based  on  (2)  is  much  simpler, 

and  is  given  in  part,  in  Appendix  B  for  locations  x  >  xc  where  xc  is 

defined  in  Figure  4a  as  the  location  of  the  liquidus  when  the  chill 

face  reaches  the  eutectic  temperature.  For  small  g  ,  condition  (2) 

E 

reduces  to  condition  (3)  which  is  by  far  the  simplest  for  analytical 
solutions.  The  equation  for  (3)  is  written: 


-  x. 


8 


XL  -  *E 


(33) 


Numerical  examples  in  Appendix  B  show  that  for  aluminum-4.5  per  cent 

v  G 

copper  alloy,  values  of  — “  calculated  from  (3)  vary  by  only  a  few 
per  cent  from  those  calculated  from  (2).  Resulting  calculated  values 
of  segregation  are  similar  (within  .02  per  cent).  Condition  (3)  is, 
therefore,  used  for  the  following. 


Figure  5:  Distribution  of  fraction  liquid  in  the  liquid-solid  region 
assuming  (a)  steady  state  solute  redistribution,  (b)  linear 
g^  versus  x'  assuming  g£  as  given  by  steady  state  equation, 

(c)  linear  gT  versus  x'  assuming  g-  ■  0. 
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2.  Macrosegregation,  x  >  x  . 

c 


Through  manipulation  of  equations  (22b),  (32),  and  (33),  an 
v  G 

expression  for  Is  readily  obtainable  (Appendix  B)  and  Is  written: 


v  G 
x 


i  (_g_ 
2^1-8 


gTno  +  2 
)(  -=-9 -  ) 

®Lnq  +  1 


where  n 


nL  ~  nE 
n„ 


(34) 


Combining  equation  (34)  with  (20)  and  Integrating  from  g^  ■  1 
at  CL  ■  CQ  to  gL,  yields  the  "solute  redistribution  expression" 
comparable  to  equations  (28),  (30)  and  (31): 


_  (1  -  a)  (1  -  k) 

-  (1  -  k)  g  n  a  +  1  a 

C.g,  t  1  (35) 


0  L 


n  a  +  1 

q 


where  a  ■  1  -  ~ 


and  macrosegregation  is  readily  determined  graphically  or  analytically 
from  (35)  and  (24  -  26). 


3.  Macrosegregation  at  x  <  x£. 

For  locations  x  <  x  ,  a  somewhat  more  complicated  expression  for 

v  G 

— — -  results  than  that  of  equation  (34).  This  expression,  derived  in 
Appendix  B  is: 


2(x  -  xg)  +  gLnq  (x  -  2x£) 
(x  -  xE)2  (qLnq  +  1) 


(36a) 


X 


E 


11  +  8L°E  (no 


St  n 
bL  q 


~  2)  +  aLgL 
+  1 


where 


34 


And  can  now  be  determined  as  a  function  of  by  combining 
equations  (36)  and  (20)  and  integrating  between  the  limits  g^  ■  1 
at  CL  ■  Cq  and  gL,  (Appendix  B) : 


1  -  k 

N  +  Mg  ~  M 
CL  "  C0  (  N  +  M  * 


(37) 


where 

Bx  (2P  -  xn  ) 

N  -  - - - 3-  (37a) 

2P 

xBOn  +  2(1  -  B)P2 

M  -  - 3 - 5 -  (37b) 

2P 


P  ■  n  x  -  a„n  +  a„  -  a. 

q  E  q  E  L 


(37c) 


Q  =  n  x  -  2a„n  +  2a.,  -  2aT  (37d) 

q  E  q  L  L 

It  should  be  noted  that  N,  M,  P  and  Q  are  constants  at  any 

distance  0  <  x'  <  x  . 

c 


Analytical  solutions  for  solute  redistribution  between  x  *  0  to 

c  c 

x  ■  x.  are  obtained  by  solving  equation  (37)  up  to  g  «  g  where  g 
L  L  L  L 

is  the  volume  fraction  liquid  at  location  x  when  x  =  0.  After 

C  h  L 

is  reached  equation  (35)  is  used  by  integrating  between  the  limits 
h  ‘  *L  8t  CL  ”  CL  and  CL-  gL 


c 


L 


(1  -  a) (1  -  k) 

g  -(1  -  k)  g  n  a  +1  a 

CL  (  -T  >  I  ~c  - ~  1  <38> 

SL  SLnqa  +  1 
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A.  Macrosegregation,  Final  Transient. 

Consider  macrosegregation  in  the  last  region  to  solidify  of  a 
unidirectional  ingot  of  length  L.  Final  composition  in  this  region 
will,  in  general,  be  less  than  CQ  (for  positive  3,  k  <  1)  as  a  result 
of  draining  away  of  solute-rich  feed  metal  to  feed  shrinkage. 

When  the  dendrite  tips  reach  the  ingot  end  (x^  ■  L) ,  at  time 
t  *  t®  subsequent  solidification  drains  liquid  back  from  the  tips, 
voids  result  since  no  further  supply  of  feed  metal  exists.  The 
location  xv,  of  the  end  of  the  liquid  at  time  t^  is: 

t 

v 

L  -  x  ■  -/  v  dt  (39) 

V  X 

t° 

V 

From  which,  by  appropriate  substitution,  equation  (39)  can  be  solved 

v 

to  determine  liquid  fraction,  at  This  liquid  fraction  will 

all  be  drained  from  the  volume  element  at  xv,  leaving  a  void  of 

y 

volume  1  -  g^»  Finally,  local  solute  redistribution  in  the  ingot  at 

v 

x^  is  the  same  for  g^  >  g^  as  in  cases  previously  considered,  but 
final  composition  is  lower  because  of  the  voids  present. 

v 

For  the  special  case  of  n^  «  0,  the  critical  value,  g^  at  x 
needed  to  calculate  macrosegregation  (from  Appendix  C) : 

L  -  x 

v 

(B)(xl  -  Xg ) 


v 


(AO) 
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5.  Numerical  Examples. 

Some  numerical  examples  of  the  foregoing  are  considered  here, 
for  unldlrectlonally  solidified  aluminum-4.5  per  cent  copper  using 
the  same  value  for  constants  as  in  previous  examples. 

(a)  Constant  Thermal  Gradient. 

Figure  6  shows  composition  versus  distance  from  the  chill  for 
constant  thermal  gradient  (n^  =  0) ,  assuming  a  total  length  of  mushy 
zone  at  steady  state  of  0.2  L  where  L  is  ingot  length.  Note  the 
initial  transient  (where  alloy  composition  is  above  C^)  is  0.2L  in 
length.  At  the  chill  face,  the  maximum  composition  calculated 
earlier  is  reached  (4.92  per  cent  copper).  The  final  transient,  result¬ 
ing  from  depletion  of  feed  metal,  is  shown  on  an  expanded  scale. 

In  later  graphs  this  zone  is  omitted  for  clarity. 

(b)  Thermal  Gradient  Decreasing  with  Time. 


Figure  7  shows  composition  versus  distance  curves  for  several 

cases  of  varying  n  .  Curve  1  shows  segregation  for  a.  and  ar  =  0, 

q  Lb 

n^  *  0.52.  (These  are  approximate  values  to  be  expected  for  aluminum- 
4.5  per  cent  copper  alloy  solidifying  against  a  metal  chill,  no 


interface  resistance  as  sketched  in  Figure  4c.)  Note  that  overall 


composition  of  the  ingot  is  imperceptibly  higher  than  CQ  (approxi¬ 
mately  4.53  as  compared  with  CQ  ■  4.50)  and  that  no  solute  rich 
"initial  transient"  zone  exists. 


r 


i 
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Figure  6:  Composition  versus  distance  from  chill  wall  for  constant 

thermal  gradient  (n_  *  0).  Length  of  mushy  zone  at  steady 
state  is  0.2L. 
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Figure  7:  Composition  versus  distance  from  chill  wall  for  varying 
thermal  conditions:  (1)  aE  ■  a-  ■  0,  nq  ■  .52, 

(2)  aL  -  0,  aE  -  -2,  nq  -  .52,  (3)  aL  «  ag  -  0,  n  -  °°, 

(4)  a,  -  0,  a£  ■*  nQ  -  .52. 
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Curve  (2)  is  identical  to  curve  1  except  that  reasonable  values 
are  employed  for  a^  and  a^  for  the  case  of  moderate  mold-metal  inter¬ 
face  resistance  to  heat  transfer  (a^  ■  0,  ag  ■  -2).  As  in  the  constant 
thermal  gradient  case  surface  composition  again  rises  to  the  maximum 
chill  face  composition  (4.92  per  cent  copper).  Composition  of  the 
bulk  ingot  is  unchanged  (4.53  per  cent  copper). 

It  is  possible  to  envlssion  two  other  simple  thermal  cases.  In 
one,  excessively  steep  thermal  gradients  are  maintained  so  n^  -►  0, 

Sjy  ■+  0,  a^  ■+  0;  here  composition  is  maintained  at  exactly  Cq  through¬ 
out.  In  a  second,  a^  and  a^  remain  zero  so  the  chill  face  is  brought 
instantaneously  to  the  eutectic  temperature  Tg.  Thereafter  the  chill 
face  is  held  at  Tg  while  the  llquidus  isotherm  moves  with  a  velocity 
proportional  to  t^  (i.e.,  n^  -*•“>).  In  this  case,  equation  (24) 
becomes : 


(41) 


and  the  solute  redistribution  equation  is  written  as: 

_  1  -  k 

CL  “  C0^8L^  3  where  a  -  1  -  (42) 


’*3crosegregation  calculated  from  the  solute  redistribution  expression 
(4 Z )  is  shown  as  curve  3,  Figure  7. 

Finally,  for  the  case  of  solidification  under  extremely  shallow 
gradients  (ar  ■  ®) ,  the  limit  of  segregation  calculated  using  equations 
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(36a)  and  (36b)  is  precisely  the  same  as  that  at  the  chill  face 

(Figure  4,  curve  4).  This  is  most  easily  seen  by  insepection  of 

equation  (20).  At  shallow  gradient,  and  finite  cooling  rate,  the 
v  G 

x 

term  — —  approaches  zero,  as  at  a  chill  face. 

(c)  Variable  Thermal  Gradients. 


Consider  unidirectional  solidification  with  movement  of  isotherms 
as  sketched  in  Figure  8.  Resulting  composition  is  shown  in  Figure  9. 
Region  A  is  the  initial  transient,  and  final  composition  is  similar 
to  that  shown  in  Figure  6.  In  Regions  B  and  F,  n^  »  0  and  final  compo¬ 
sition  is  Cq.  Region  D  is  drawn  with  n^  =  2.5  with  the  resulting 
increase  in  composition  (see  Figure  10) .  Regions  C  and  E  are  transient 
zones  between  B-D  and  D-F  respectively.  The  compositional  variations 
in  these  two  transient  regions  are  depicted  schematically. 


Macrosegregation  in  regions  B,  D,  and  F  is  readily  calculated 
from  the  solute  redistribution  equation  (35)  shown  as  Figure  10. 

Figure  9  shows  the  final  macrosegregation  in  the  ingot.  Note  the 
positive  segregation  (solute  rich  "band*  )  at  region  C  resulting  from 
decreasing  thermal  gradient.  Similarly  an  increasing  thermal  gradient 
will  result  in  negative  segregation  (solute  poor  "band"). 

(d)  Variable  Ingot  Cross  Section. 


Consider  an  ingot  of  design  shown  ii.  Figure  11  solidifying 
unidirectionally  upwards  with  constant  thermal  gradient  for  liquid- 
solid  region  is  of  length  0.2L  for  x^.  >  0.  The  ingot  is  reduced, 
on  an  area  basis  10/1  ai.  0.5L. 
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Figure  8:  Isotherm  movement  for  a  unidirectionally  solidified  ingot. 
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Figure  9: 


Composition  versus  distance  from  the  chill  wall  for  Ingot 
with  isotherm  movement  as  shown  in  Figure  8. 
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Final  solid  composition,  Figure  12,  Is  as  in  the  comparable 

constant  cross  section  Ingot  (Figure  6)  up  to  x/L  -0.5.  At  x/L  - 

v  G 

0.5,  v  ,  and  hence  Increases  discontlnuously  by  a  lactor  of  10 

(see  equation  23).  Hence,  solid  composition  at  x/L  ■  0.5  is  given 

VxG  1  -  B 

directly  by  Figure  3  for  — —  — - —  -  10  (Cg  -  2.3  per  cent  copper). 


For  the  transient  zone,  0.5  <  x/L  <  0.7  a  region  exists  where 

locax  solute  redistribution  is  expressed,  from  equation  (20),  for 

c  .  V::G  1-6  in  . 

gL  >  gL  where  — - —  -  10  by: 


g 


L 


1  +  9B 
1  -  k 


(43a) 


V  ,  q 

For  g.  <  gf,  C.  >  cf,  — — - —  -  1  and  equation  (20)  in  the 

L  L  L  L  G  D 

integrated  form  becomes 


(43b) 


For  0.7  <  x/L  <  1.0  final  composition  in  the  reduced  cross 

section  ingot,  Figure  12,  is  again  Identical  to  that  of  the  comparable 

v  G 

constant  cross  section  case.  Figure  6,  where  (1  -  6)  -  1. 


A  second  type  of  ingot  of  variable  cross  section  is  one  in  which 
cross  section  varies  continuously.  Consider,  for  example,  the  prism 
of  Figure  13  in  which 


A' 


(44) 
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Figure  .12:  Composition  versus  distance  from  chill  for  Ingot  shown  in 
Figure  11. 
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where:  A'  ■  cross  sectional  area  at  x' 

Aq  ■  cross  section  area  at  chill  face 

v  G 

Again  assuming  constant  temperature  gradient,  —j—  is  calculated 
from  equation  (23)  which,  for  the  assumptions  used  herein  reduces 
to,  for  x  >  x£  (Appendix  D) 


V  .  b  -  l  ..  2U*l  -  V  *  2xar+ 

e  1-6'L-x'  2(x  -  x_) 

L  E 


(45) 


where  aT  »  a.  -  a„ 

T  L  E 

Local  solute  distiibution  and  macrosegregation  are  then  calculated 
as  before  by  integrating  (20)  after  substituting  (45)  and  then 
solving  ( 26) . 

For  x  <  x£  equation  (45)  becomes 


Vx^  8  .  1  v  ,  x(2L  -  x)  , 

e  1  -  8  (  L  -  x  )[  '  2gL  J 


(45b) 


and  macrosegregation,  is  determined  by  a  procedure  similar  to  that  for 
the  ingot  of  constant  cross  section  (equations  36  -  38). 

For  g^  >  g£,  <  c£  equation  (45b)  is  substituted  in  (20)  and 

integrated  to  solve  for  solute  redistribution.  When  g  <  g^,  C  >  CC 

L  L  L  L 

equation  (45a)  is  used  in  (20)  and  integrated  between  appropriate 
limits.  The  results  are  plotted  as  Figure  14. 
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(e)  Variable  k,  6. 

To  permit  analytic  solutions,  the  foregoing  calculations  have 
been  made  assuming  constant  k  and  constant  8.  In  aluminum-4.5  per 
cent  copper  these  factors  are,  however,  not  constant.  Figure  15  shows 
k  versus  C^;  Figure  16  summarizes  best  available  data  on  Pg  and  p^ 
versus  C^;  Figure  17  shows  6  versus  C^»  calculated  from  the  data  of 
Figure  15.  Note  that  in  this  alloy,  $  is  negative  during  much  of 
solidification  (negative  6  denotes  expansion  on  solidification). 
However,  a  net  flow  of  feed  metal,  negative  v  ,  still  results  through- 

X 

out  steady  state  solidification,  because  of  thermal  contraction  of 
the  liquid  (i.e.,  the  last  term  on  the  right  of  equation  (12)  is 
sufficiently  large  (negative)  that  v^  is  negative  even  though  the 
first  term  on  the  right  is  positive). 

Figure  18  shows  local  solute  redistribution  at  the  chill  face 
for  aluminum-4.5  per  cent  copper,  determined  using  the  data  of  Figures 
15  -  17  by  graphical  integration  of  equation  (20).  For  comparison 
with  the  case  of  constant  8  (and  hence  constant  Pg,  p^) ,  the  data  are 
plotted  as  weight  fraction  solid,  fg,  where: 

fs  “  PS(1  -  «£>  ‘  (4M 

where:  Pg  and  p£  ■  densities  of  the  a  phase  and  eutectic 

respectively  (at  just  below  T^) 

gg  ■  volume  fraction  eutectic 


0.20 
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Figure  15:  Equilibrium  partition  ratio  versus  liquid  composition  for  aluminum  rich,  aluminum- 
copper  alloy. 
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WEIGHT  FRACTION  SOLID  (fs) 


Figure  18:  Steady  state  solute  redistribution  for  two  cases: 

(a)  constant,  k,  8  direct  calculation  of  weight  -  volume 
fraction  solid,  and  (b)  variable  k,  8  -  calculation  of 
volume  fraction,  from  density  data  Figure  16  conversion 
to  weight  fraction. 
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Graphical  integration  of  curve  (b) ,  Fieure  18,  gives  a  final 
chill  face  composition  of  5.03  per  cent.  This  curve,  based  on  best 
available  data  for  alloy  constants,  yields  a  final  value  of  chill 
face  composition  close  to  that  calculated  by  Schell ^  using  slightly 
different  constants.  The  curve  (a)  in  Figure  18  was  calculated 
with  constant  k  (k  ■  .172),  and  constant  B  (8  ■  .055),  with  the 
value  of  8  chosen  so  that  final  chill  surface  composition  is  that 
calculated  by  Schell  (4.93).  Note  the  similarity  of  curves  (a)  and 
(b) ,  making  reasonable  the  approximation  used  throughout  much  of 
this  work,  of  constant  B. 

PART  III:  NON-UNIDIRECTIONAL  HEAT  FLOW 

A.  Introduction 

The  foregoing  has  been  concerned  only  with  cases  where  isotherms 
could  be  assumed  planar,  and  heat  flow  unidirectional  from  the  metal 
to  a  single  flat  mold  wall.  In  this  section  we  examine  macrosegrega¬ 
tion  in  castings  or  ingots  in  which  heat  flow  is  not  unidirectional. 
Types  of  segregation  include  (1)  centerline  segregation,  (2)  negative 
"cone"  of  segregation,  (3)  positive  segregation  beneath  risers  (hot 
tops),  and  (4)  others.  It  is  assumed  in  most  of  the  following,  as 
before,  that  feeding  is  complete  (no  voids  form  during  solidifica¬ 
tion).  This  however,  is  expected  to  be  an  increasingly  poor 
assumption  as  thermal  gradients  decrease. 


Figure  19  shows  schematically  solidification  of  an  alloy  in  a 
metal  mold.  In  Figure  19a,  a  fully  liquid  region  remains  at  the  ingot 
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center  and  in  Figure  19b,  at  a  later  stage  of  solidification,  this 
fully  liquid  zone  has  disappeared.  It  will  be  shown  below  that  some 
centerline  segregation  will,  in  general,  result  from  flow  to  feed 
shrinkage  during  and  after  the  time  of  solidification  sketched  In 
Figure  19b.  Also  it  will  be  shown  that  flow  to  feed  shrinkage 
throughout  most  of  solidification  must  result  in  negative  segregation 
near  the  ingot  base  (negative  "cone"  of  segregation) . 

Figure  20  is  a  sketch  of  a  plate-type  sand  casting  which  may  be 
made  with,  or  without,  a  metal  chill  at  the  end  of  the  plate  opposite 
the  riser.  Figure  21  shows  schematically  the  solidification  of 
aluminum-4.5  per  cent  copper  alloy  in  a  3/4- inch  thick  plate  of  this 
geometry.  Temperature  differences  throughout  most  of  the  plate 
during  solidification  are  small  compared  with  the  freezing  range  of 
the  alloy,  and  so  solidification  proceeds  in  "mushy"  manner  as  shown 
in  Figure  21a.  When  a  chill  is  placed  at  the  end  opposite  the  riser, 
longitudinal  thermal  gradients  are  introduced  in  the  metal  and 
solidification  proceeds  directionally  towards  the  riser  as  sketched 
in  Figure  21b.  The  two  different  types  of  solidification  result  in 
two  distinctly  different  patterns  of  macrosegregation  as  discussed 
below. 

As  a  final  example  of  macrosegregation  in  non-unidirectional 
solidification,  the  case  of  segregation  under  a  riser  ("hot  top")  of 
insufficient  size  is  considered.  As  example  we  will  examine  a  cube, 
rir.ered  with  a  tall  riser  of  square  cross  section,  which  is  smaller 
than  that  of  the  cube.  Figure  22.  The  small  riser  cools  through  the 
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Figure  19a:  Sketch  of  ingot  solidification  In  a  metal  mold 
showing  liquid,  liquid  plus  solid,  and  fully 
solid  regions. 
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Figure  19b:  Sketch  of  Ingot  solidification  in  a  metal  mold 
during  late  stages  of  solidification,  after 
llquldus  isotherms  have  met  at  the  casting 
centerline. 
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solidification  range  somewhat  more  rapidly,  than  the  casting  and 
therefore  contains  somewhat  more  solid  during  most  of  solidifica¬ 
tion  than  does  the  casting.  It  will  be  shown  below  that  this 
results  in  positive  macrosegregation  under  the  riser. 

B.  Centerline  Segregation 

1.  No  Contraction  Strain. 

Consider  solidification  from  two  flat  parallel  faces,  so 

isotherms  are  planar  and  meet  at  the  centerline,  Figure  19.  Before 

the  centerline  reaches  the  liquidus  temperature,  as  in  Figure  19a, 

liquid  flow  is  expected  to  occur  in  the  y  vertical  direction  only  in 

the  region  of  no  solid  (region  between  the  dendrite  tips  and  the 

centerline)  within  the  liquid-solid  region,  v  ■  v  "0,  and  v  is 

y  z  x 

as  in  steady  state  unidirectional  solidification. 

When  the  liquidus  isotherm  reaches  the  centerline,  Figure  19b,  flow 

within  the  liquid-solid  region  must  occur  in  the  y  direction.  Three 

possible  flow  paths,  with  resulting  segregation,  will  now  be  considered. 

In  each,  for  simplicity,  assumptions  employed  will  include  constant 

B,  k,  G,  c,  and  linear  distribution  of  solid  in  the  liquid-solid  zone. 

Alloy  used  as  example  is  again  aluminum-4.5  per  cent  copper;  n  is 

9 

assumed  equal  to  0  and  at  steady  state  the  thickness  of  the  mushy 
zone  is  0.2L,  where  L  is  half  thickness  of  this  ingot. 

(a)  Vy  *  0  within  the  liquid-solid  zone  except  at  x'  ■  L.  Here, 
all  flow  in  the  y  direction  is  in  a  layer  of  vanishingly  small  thickness 
at  the  centerline.  Hence,  v^  at  all  locations  up  to  x'  ■  L  is  exactly 
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as  in  steady  state  solidification,  and  G/t  is  constant.  Final  local 
solute  content  is  as  in  steady  state  solidification  and  no  macrosegre- 
gation  results  at  the  centerline,  curve  (1),  Figure  23. 


(b)  v  ■  0  at  x  <  L;  v  “0  after  x  »  L.  Here,  the  instant 
y  LX  L 

the  llquidus  isotherms  reach  the  centerline,  flow  in  the  x  direction 

is  assumed  to  cease.  If  at  this  instant  the  fraction  liquid  at  x'  is 

c  c 

g^»  then  it  is  clear  that  for  >  g^,  solute  redistribution  is  as  in 

steady  state  solidification,  and  for  g^  <  g^,  solute  redistribution 

is  at  a  chill  face.  Hence,  local  solute  redistribution  is  calculated 

from: 


8t 


C.  1  -  k 

<r> 

L0 


f°r  gL  >  gL 


(45a) 


8t 


c  .  JL  . 


1  -  B 
1  -  k 


for  8L  <  8L 


(45b) 


and  segregation  is  calculated  as  before,  and  shown  as  curve  (2), 
Figure  23. 


(c)  Vy  *  0  at  x^  <  L;  v^  ■  f(g^)  ^or  xl  >  cases  (1) 

and  (2)  represent  limits  on  centerline  segregation  which  can  occur 
in  the  examples  chosen,  v  ,  and  hence  v  and  segregation  can  generally 

y  X 

be  expected  to  lie  between  the  extremes  chosen.  For  example  when 
Vy  *  f ( gj^)  solute  content  in  the  final  transient  region  will  at  all 
locations  lie  between  case  (1)  and  case  (2).  At  the  beginning  of  the 


transient  region  composition  will  be  while  at  the  centerline  it  is 


the  same  composition  as  case  (2)  where  no  flow  has  occurred  in  the  x 
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(b) 


Figure  23:  Final  composition  versus  distance  for  two  cases  of  Ingot 
mold  solidification,  (a)  Steady  state  Vy  ■  0,  within 
the  liquid-solid  zone  except  at  x'  ■  L.  (b)  Steady  state 
vy  -  0  at  xL  <  L;  vx  -  0  after  xL  -  L. 
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direction.  The  shape  of  the  curve  in  between  these  two  points, 
however,  is  determined  by  the  dependence  of  v^  on  gj  .  Solution 
requires  detailed  consideration  of  fluid  flow  within  the  porous 
medium. 


2.  Contraction  Strain. 

In  the  foregoing,  movement  of  solid  during  solidification  has  not 
been  considered.  Such  movement  might  occur  by  (a)  floatation  or 
settling,  (b)  contraction  or  other  strain.  Contraction  strains  can 
result  in  internal  fissuring  and/or  "filled  hot  tears"  at  any  loca¬ 
tion,  x,  y,  z  in  an  ingot;  one  example  of  this  is  at  the  ingot 
centerline,  considered  below. 

Consider  the  case  examined  in  the  preceding  section  of  solidifica¬ 
tion  with  planar  isotherms  from  two  flat  parallel  faces  meeting  at  a 
centerline  as  shown  in  Figure  19 r  except  that  "fissuring"  will  occur 
from  thermal  contraction  of  the  already  solidified  metal,  at  the 
casting  centerline,  x'  ■  L.  When  the  centerline  just  reaches  the 
liquidus  temperature,  solid  begins  to  form  and  is  drawn  away  from  the 
centerline  at  a  velocity  v  from  thermal  contraction;  hence,  the 
"volume  element"  associated  with  this  solid  is  also  drawn  away  at 
velocity  v'.  At  subsequent  temperatures  T  below  the  liquidus,  other 
volume  elements  form  at  the  centerline,  in  liquid  of  composition 
C“  >  CQ  (for  k  <  1). 

For  the  volume  element  which  forms  at  the  centerline  at  temperature 
T  below  the  liquidus  (when  the  liquid  composition  at  the  centerline  is 


66 


raised  to  C®) ,  equation  (20)  Is  written: 


.  -  -  t  ( 


a  vvG  dCL 

T><  1+T-> 


The  final  position  of  the  volume  element  at  the  end  of 
solidification  (t^  +  t^)  is: 


*L  +  tf 


f  v'dt 


where:  t^  and  t^  are  times  for  initiation  and  completion  of 

solidification  at  x 

A^  ■  distance  from  casting  centerline 


If  the  only  strain  is  from  thermal  contraction  and  the  strain  is 


small: 


x  *  x 


/  Tdx’ 


x  =  0 


where:  T  *  temperature  at  x 

a  *  linear  coefficient  of  thermal  contraction  at  T 


Qualitative  conclusions  from  the  foregoing  are  that  thermal 
contraction  can  cause  segregation  at  the  centerline  (or  elsewhere) 
in  the  solidifying  ingot.  Segregation  resulting  from  this  contraction 


is  additive  to  that  resulting  from  solidification  contraction. 
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Equations  (46)  and  (48)  can  be  used  to  evaluate  thickness  of 
segregated  zone,  and  g^  and  at  any  location  within  the  zone  during 
solidification.  Hence,  final  segregation  at  the  centerline  is  readily 
calculated  as  above. 


Consider,  for  example,  the  unidirectionally  solidification  aluminum- 
4.5  per  cent  copper  alloy  used  as  example  above,  constant  k,  6,  G,  e 
linear  distribution  of  solid  in  the  liquid-solid  zone  (with  z) . 

Assume  near  the  region  in  question  fluid  flow  (to  feed  thermal  contrac¬ 
tion  and  shrinkage)  occurs  only  by  (a)  flow  down  the  ingot  centerline, 
parallel  to  isotherms,  and  (b)  flow  perpendicular  to  isotherms  at 
locations  other  than  the  centerline. 


Hence, 


v  G 
x 


8 


1  -  8 


so  equation  (31)  becomes: 


r°  -  —  ■  — 

(It  )  k 

CL 


(49) 


Figure  24  plots  Cg  versus  gg  for  three  values  of  C®,  (a)  C®  ■  Cq  = 
4.5  per  cent  copper,  (b)  C£  *  20  per  cent  copper,  and  (c)  C®  ■  33  per 
cent  copper. 

Assume  now  thermal  contraction  was  such  that  each  volume  element 
moved  .001  cm/per  cent  increase  in  solute  at  the  centerline.  Hence, 
a  segregated  centerline  results  of  thickness  (.001) (33)  ■  .033  cm  in 
which  the  average  composition  varies  linearly  from  4.5  per  cent  copper 
to  33  per  cent,  Figure  25. 


COMPOSITION  SOLID  l PERCENT  COPPER) 
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VOLUME  FRACTION  SOLID  (gs) 


Figure  24:  Solid  composition,  C-,  versus  volume  fraction  solid,  g_, 
for  three  values  of  C£;  (a)  C£  ■  Cq  ■  4.5  per  cent 
copper,  (b)  C£  ■  20  per  cent  copper,  (c)  C£  ■  33  per 
cent  copper. 
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4  3  2  1 

DISTANCE XIO3  (CM.) 


? 


CENTERLINE 


Figure  25:  Centerline  segregation  resulting  from  thermal 

con' raction  of  .001  cm/per  cent  increase  in  solute 
at  ‘.he  centerline. 
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C.  Negative  "Cone11  of  Segregation  at  Ingot  Base 

Consider  solidification  of  the  bottom  portion  of  the  ingot 
illustrated  in  Figure  19.  Assume  as  limiting  case  that 

(1)  at  the  periods  of  solidification  represented  by  Figure  19a, 
feed  metai  flows  into  the  bottom  "cone"  (represented  by  the 
dendrites  growing  vertically)  only  from  the  fully  liquid 
region,  and 

(2)  after  the  dendrites  growing  from  the  sidewalls  reach  the 
centerline,  flow  into  the  "cone"  is  all  through  a  small 
channel  at  the  apex  of  the  cone. 

Now,  assume  an  ingot  of  long  rectangular  cross  section  (long 
dimension  perpendicular  the  planr  of  the  paper).  The  "cone"  then 
is  a  prism,  and  hereafter  referred  to  as  such.  Moreover,  fluid 
flow  into  this  prism  is  exactly  as  in  the  prism  ingot 
discussed  earlier  (Figure  13)  and  final  composition  within  the 
prism  is  expected  to  be  as  given  in  Figure  14.  Note  that  composi¬ 
tions  particularly  near  the  apex,  are  substantially  below  nominal, 
in  qualitative  agreement  with  experiments  on  large  ingots  which  show 
what  is  commonly  termed  the  "negative  cone"  of  segregation. 

D.  Macrosegregation  in  Sand  Cast  Plates 

In  sand  castings  temperature  gradients  within  the  solidifying 

metal  are  generally  small,  because  of  the  low  heat  diffusivity  of 

the  mold  compared  with  that  of  the  metal.  Fhere  these  are  negligible, 

v  G 

and  where  feeding  is  complete,  the  term  in  equation  (20)  is  zero. 
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For  the  plate  shown  in  Figure  20  solidification  is  sketched  in 
Figure  21.  Significant  macrosegregation  is  expected  through  the 
plate,  approaching  the  maximum  expected  at  a  chill  face  in  unidirec¬ 
tional  solidification  (as  temperature  gradient  in  the  longitudinal 
direction  goes  to  zero,  assuming  complete  feeding).  In  the  region 
near  the  riser,  some  thermal  gradients  are  expected,  and  hence,  some 
reduction  in  final  composition  in  this  area  is  also  expected, 

Figure  26  illustrates  this  schematically  for  aluminum-4.5  per  cent 
copper  alloy. 

In  the  end  chilled  plate.  Figure  21b,  steep  longitudinal  gradients 
are  Introduced  by  the  chill,  and  although  some  heat  is  lost  to  the 
cold  sand  mold,  isotherms  are  expected  to  be  essentially  planar,  and 
perpendicular  to  the  longitudinal  (x')  direction.  Since  the  isotherms 
are  assumed  planar,  relations  developed  in  earlier  sections  for 
"unidirectional  solidification"  apply  here  as  well.  Movement  of  the 
isotherms,  however,  is  expected  to  depend  on  plate  geometry  and  heat 
diffusivity  of  the  mold  as  well  as  on  rate  of  heat  removal  through 
the  chill.  Figure  26  illustrates  composition  distribution  expected 
assuming  (1)  movement  of  isotherms  throughout  most  of  the  plate  is  as 
in  unidirectional  solidification  from  a  metal  chill,  no  interface 
resistance,  but  (2)  isotherms  approach  each  other  near  the  riser  end 


of  the  plate. 
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E.  Under  Riser  Segregation  Positive  Segregation  Near  "Hot  Top" 

For  one  example  of  under  riser  segregation  consider  the  risered 
sand  cast  cube  of  Figure  22.  Here,  the  riser  chosen  is  inadequate 
to  feed  the  cube  and  at  any  time  during  solidification  the  metal  in 
the  riser  area  will  contain  more  solid  than  the  metal  in  the  casting, 
Figure  27. 

Thermal  gradients  are  assumed  small  in  the  horizontal  (y  and  z) 
directions  of  casting  and  riser.  The  thermal  gradient  is  also 
assumed  small  in  the  x  direction  except  in  the  region  between  the 
casting  and  riser.  To  illustrate  the  problem  simply,  assume 
temperature  difference  between  bulk  riser  and  bulk  casting  are 
constant  duiing  solidification,  and  the  temperature  gradient 
in  the  vicinity  of  the  casting  and  riser  junction  is  constant. 

Consider  first  solidification  ir>  the  bulk  riser,  well  away  from 

the  casting  riser  junction;  here,  G  «  0.  Since  v  and  e  are  finite 
v  G 

x 

— “  *  0  and  maximum  segregation,  as  at  a  chill  face,  is  expected 
(provided  no  pores  form) ,  except  in  the  region  at  the  top  of  the 
riser  where  negative  segregation  occurs  as  a  result  of  liquid  depletion. 

Similarly,  while  feeding  occurs,  solute  redistribution  in  the 
bulk  casting  must  be  as  at  a  chill  face  (at  locations  well  away  from 
the  casting-riser  junction).  However,  at  some  critical  fraction 
solid,  gg  which  must  occur  not  later  than  when  the  riser  reaches  T^, 
its  non-equilibrium  solidus,  all  feeding  is  expected  to  cease.  Final 
composition  in  the  bulk  casting  is  therefore  expected  to  be  less  than 
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Figure  27:  Sketch  of  solidification  for  cube  casting  of  Figure  22. 
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the  chill  face  maximum,  by  an  amount  determined  by  the  time  during 
solidification  when  feeding  ceases.  This  composition  will  however, 
be  at  least  slightly  greater  than  CQ. 

In  the  region  at,  or  just  below,  the  riser-casting  junction,  we 

expect  composition  to  be  significantly  higher  than  in  the  bulk 

casting.  This  ii  because,  as  long  as  flow  occurs  in  this  region,  it 

v  G 

X 

is  flowing  up  a  temperature  gradient;  hence  — —  is  positive  and 
compositions  in  excess  of  the  chill  face  maximum  are  possible. 
Compositions  in  excess  of  the  chill  face  maximum  must  result  as  long 
as  there  is  a  net  flow  of  solute  up  a  temperature  gradient,  and  this 
flow  continues  to  the  last  stages  of  solidification. 


■Wl  1'  IP.JJU. 
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III.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

In  this  section,  results  are  presented  from  experiments  designed 

to  illustrate  and  substantiate  certain  aspects  of  the  general 

macrosegregation  expression  previously  developed.  In  particular, 

v  G 

consideration  is  given  to  the  parameter,  — —  in  equation  (20) 


v  G 

)<  1  +  -r> 


(20) 


This  solute  redistribution  expression  is  for  the  case  of 
unidirectional  heat  flow,  but  where  fluid  flow  is  not  necessarily 
unidirectional. 

As  previously  discussed  for  a  particular  alloy  system  and 

composition,  values  can  be  assigned  which  effectively  average  k  and 

B  over  the  complete  solidification  range,  even  though  they  are 

usually  functions  of  C^.  Using  the  values  B  *  .055  and  k  =  .172  as 

the  constants  for  this  particular  alloy  composition,  consideration 

can  then  be  given  to  the  effect  of  solidification  variables  on  the 
v  G 

x 

parameter  — ““  .  The  following  experimental  results  are  presented 
to  show  how  variations  in  this  parameter  can  cause  compositional 
changes  throughout  a  casting  or  ingot. 

A.  Unidirectional  Heat  and  Fluid  Flow 
1.  Constant  Cross  Section. 

Consideration  is  first  given  to  experimental  results  on  four 
ingots  of  constant  cross  section  where  both  heat  and  fluid  flow  are 
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unidirectional.  The  ingots  were  made  using  the  experimental  set  up 
described  in  Appendix  E.  Ingots  1-3  are  of  square  cross  section 
7.5  cm  per  side,  20  cm  high.  Ingot  number  4  is  circular  11.5  cm  in 
diameter  and  20  cm  high. 

(a)  Ingot  Number  1. 

Ingot  number  1  was  cast  and  solidified  on  a  water  cooled  chill. 
Furnace  temperature  was  maintained  high  enough  to  obtain  a  completely 
columnar  casting  structure.  Cooling  curves  for  the  ingot  are  shown 
in  Figure  28.  From  these  curves  Figure  29  was  constructed,  showing 
liquidus  and  eutectic  isotherm  movement  as  «  function  of  time.  It 
is  interesting  to  note  that  even  though  a  water  cooled  chill  was 
used,  after  720  seconds  the  mushy  zone  extends  the  entire  length  of 
the  casting.  That  is,  at  the  same  time  the  metal  at  the  chill  first 
becomes  completely  solid,  initial  solid  is  forming  at  the  top.  This 
means  that  at  least  the  initial  stages  of  solidification,  at  all  loca¬ 
tions,  is  under  the  initial  transient  conditions  previously 
considered  (i.e.,  x  <  x£) . 

Calculation  of  composition  up  the  length  of  the  ingot  using  only 
thermal  data  was  done  as  described  in  Appendix  G.  Agreement  between 
this  and  the  experimentally  determined  curve  (using  fluorescent 
analysis  for  chemical  analysis,  see  Appendix  F)  is  shown  in  Figure  30. 

Calculated  agreement  becomes  poor  at  distances  greater  than  10 
centimeters  from  the  chill.  It  can  be  seen  from  the  thermal  data  of 
Figure  28  that  at  this  location  gradients  become  shallow  (typically 
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Figure  29:  Liquidus  and  eutectic  isotherm  movement  for  Ingot  No. 
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Figure  30:  Experimental  and  calculated  (from  thermal  data)  solid  composition  versus  distance  from 
chill  in  Ingot  No.  1. 
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2  -  5°C/cm),  a  condition  which  can  lead  to  porosity  and  therefore 
lower  composition. 

(b)  Ingot  Number  2. 

Ingot  number  2  was  also  cast  and  solidified  against  a  water 
cooled  chill.  Furnace  temperature  was  lower  than  for  the  preceding 
ingot,  and  resulted  in  an  equiaxed  ingot  structure,  even  though  heat 
flow  was  unidirectional.  Insufficient  thermal  data  is  available  for 
complete  analysis,  but  liquldus  and  solidus  isotherm  movements  were 
recorded.  An  initial  transient  of  about  5  centimeters  (Figure  31) 
exists  which  is  seen  from  the  increased  composition  at  locations 
less  than  5  cm.  Figure  32. 

In  the  region  after  the  initial  transient,  n^  (equation  34), 
which  reflects  the  relationship  between  isotherms,  changes  very  little. 
At  locations  greater  than  Z  cm,  n^  varies  from  .8  to  1.0.  Figure  10 
shows  that  within  this  range  of  n^  values  composition  is  not  expected 
to  vary  more  than  .03  per  cent  and  within  experimental  limits  should 
be  relatively  constant.  This  is  the  case  in  the  region  from  5  to 
about  15  cm  from  the  chill  where  composition  is  about  4.68  per  cent 
copper.  It  is  also  expected  (Figure  10)  that  within  this  region 
(i.e.,  5-15  cm,  n  ■  .8  -  1.0)  average  composition  will  be  slightly 
greater  than  CQ.  For  this  ingot  is  4.60  per  cent  copper,  less  than 
the  composition  in  the  5  -  15  cm  range  (4.68  per  cent).  Toward  the 
top  of  casting  (greater  than  15  cm)  composition  begins  to  drop  off 


since  solute  must  be  conserved. 
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Figure  31:  Liquidus  and  eutectic  isotherm  movement  for  Ingo 
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(c)  Ingot  Number  3. 

Ingot  number  3  was  cast  and  solidified  against  a  chill  which  had 
both  air  and  water  for  a  coolant.  Figure  33  presents  the  experimental 
and  calculated  (from  thermal  data)  composition  curves.  At  locations 
greater  than  15  centimeters  where  gradients  are  shallow  and  porosity 
forms,  agreement  between  experimental  and  calculated  curves  breaks 
down.  The  general  method  used  for  calculating  the  curve  is  outlined 
in  Appendix  G. 

(d)  Ingot  Number  4. 

Ingot  number  4  was  cast  against  a  chill  and  air  cooled  to  start 
solidification.  After  the  liquidus  had  extended  2.2  -  2.5  cm  into 
the  casting,  before  the  liquid  at  the  chill  reached  the  eutectic 
temperature,  the  air  was  replaced  by  a  water  coolant.  The  resulting 
thermal  curves  for  some  typical  isotherms  are  shown  in  Figure  34. 
Resulting  composition  near  the  chill  face  is  shown  in  Figure  35. 

When  comparison  is  made  between  the  liquidus  and  eutectic  isotherms 

the  initial  transient  (x  <  x  )  extends  almost  4  cm  into  the  casting. 

c 

This  precludes  quantitative  comparison  in  this  region  of  the  liquidus 
and  eutectic  isotherm  slopes  to  obtain  n^  because  use  of  equation  (34) 
requires  x  >  x^.  Also,  complete  analysis  as  in  ingots  1  and  3  was  not 
possible  due  to  insufficient  thermal  data  in  the  critical  region. 
However  a  qualitative  comparison  can  be  made  between  the  liquidus 
isotherm  and  an  isotherm  less  than  100  per  cent  aolid  (eutectic) .  For 
example  consider  the  liquidus  and  30  per  cent  solid  isotherms  in 
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DISTANCE  FROM  CHILL- WALL  (CM.) 

Figure  35:  Solid  composition  versus  distance  from  the  chill  in  Ingot  No. 
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Figure  34.  We  can  define  n  as: 

qx 

nT  -  n 
L  x 

n  ■  •  ~~ 

qx 

where  n  «  slope  of  an  isotherm  movement  line  where 
0  <  x  <  100  per  cent 

Now  n  can  be  calculated  which  will  give  a  qualitative  indication 
of  whether  compositions  should  be  greater  or  less  than  Cq  even  though 
quantitative  results  are  not  available.  In  the  region  where  x  is  less 
than  2.5  cm  comparison  between  the  liquidus  and  the  30  per  cent  solid 
isotherm  slopes  gives  n^x  as  -.9.  The  maximum  negative  deviation 
possible  between  two  isotherms  is  -1,  therefore  significant  negative 
departure  from  is  expected  in  the  region  up  to  2.5  cm.  After 
x  ■  2.5  cm  the  relationship  between  the  liquidus  and  all  other  slopes 
becomes  about  the  same;  that  is,  n^  becomes  constant  which  predicts 
a  constant  composition  (Figure  35). 

2.  Reduced  Cross  Section. 

The  next  group  of  castings  were  reduced  in  cross  section  as  shown 
in  the  schematic  diagram,  Figure  11.  While  the  castings  will  be 
considered  again  in  a  later  section,  attention  in  this  portion  is 
directed  only  at  the  effect  of  the  reduction  on  composition  through 
the  central  portion  of  the  ingot.  That  is  consideration  is  given  to 
a  section  3/4  inch  by  3/4  inch  extending  from  the  chill  surface  to 
the  top  of  the  casting  through  the  geometric  center. 
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Figure  36  shows  the  composition  variations  up  the  length  of 
an  ingot  in  the  center  section  where  area  reduction  is  9/1.  As 
expected  there  is  a  sharp  minimum  at  the  intersection  of  the  large 
and  small  sections  (see  theoretical  section,  Figure  12,  variable 
ingot  cross  section).  To  the  left  of  the  intersection  there  is  no 
change  from  00  to  zero  slope  as  in  the  theoretical  case  (Figure  12). 
This  is  because  the  theoretical  model  assumes  constant  v^  at  any 
location  y  from  the  chill.  In  the  actual  ingot,  just  below  the 
reduction,  v^  will  be  greatest  at  the  center  of  the  ingot  (i.e. 
absolute  value)  and  decrease  toward  the  edges.  This  will  be 
considered  further  in  a  later  ingot,  where  extensive  thermal  data  was 
available. 

An  experiment  was  next  conducted  similar  to  the  previous,  but 
with  the  area  change  2/1;  results  are  presented  in  Figure  37.  For 
this  ingot,  P  equals  2  (Figure  3)  with  a  predicted  minimum  composition 
at  the  intersection  about  4.1  per  cent  copper.  Experimental  value  is 
4.16  per  cent.  The  general  shape  of  the  composition  curve  is  similar 
to  that  of  the  previous  figure,  except  for  a  secondary  minimum  about 
12  cm.  This  may  be  related  to  a  sudden  velocity  increase  at  this 
location  when  the  final  liquid  solidifies  just  below  the  intersection. 
This,  however,  would  require  verification  by  thermal  analysis  and  has 
not  been  done. 

The  final  ingot  presented  in  this  section  is  a  repeat  of  the  9/1 
reduction  ingot  made  with  thermocouples  cast  in  place  so  that 
experimental  and  theoretical  curves  could  be  compared  with  curves 
calculated  from  the  thermal  data.  Figure  38  shows  the  measured 
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Figure  37:  Solid  composition  versus  distance  from  the  chill  in  an  ingot  of 
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composition  using  fluorescent  analysis •  Figure  39  shows  the 
theoretical  curve  for  a  A, 5  per  cent  alloy  solidifying  with  a  mushy 
zone  0.2L  at  steady  state*  This  curve  had  previously  been  presented 
in  Figure  12.  The  curve  as  presented  in  Figure  39b  was  calculated 
using  thermal  data  from  the  actual  ingot,  as  described  in  Appendix  G. 
The  sharp  breaks  at  the  reduction  in  area  of  the  theoretical  and 
calculated  curves  are  due  to  the  assumption  that  v^  is  constant  at 
a  given  time  over  the  entire  y-z  plane  at  a  given  distance  x  from 
the  chill.  This  requires  discontinuous  change  in  v  at  the  reduction. 

X 

In  the  actual  case  this  is  not  to  be  expected  since  it  would  require, 

just  below  the  reduction,  that  feed  metal  flow  in  a  thin  layer  across 

the  top  of  the  large  section  and  then  down,  so  that  v  and  v  equal 

y  z 

zero  at  all  locations  except  precisely  at  the  junction.  Experimental 
data  from  this  ingot,  presented  in  Figure  AO,  show  that  final  solid 
composition  just  below  the  junction  is  much  lower  near  the  center  of 
the  casting  than  it  is  toward  the  edges.  This  clearly  results  from 
finite  v  and  v  ,  and  is  discussed  in  a  later  section. 

y  z 

As  a  final  example  of  the  effect  of  section  reduction  on 
composition,  partial  results  of  an  ingot  with  metal  plates  to  divert 
fluid  flow  are  presented  in  Figure  Al.  The  plates  are  spaced  so  they 
extend  half  the  width  into  the  casting  (total  width  7.5  cm).  The 
side  of  the  casting  with  the  plates  will  be  considered  in  a  following 
section.  On  the  slue  with  no  plates,  the  minimum  composition  occurs 
at  approxir ately  k  cm  from  the  chill.  The  2/1  reduction  in  area  from 
the  inserted  plates  occurs  at  5  cm.  That  the  expected  minimum  occurs 
below  the  reduction  seems  possible  since  from  fluid  flow  considerations 
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Figure  38:  Solid  composition  versus  distance  from  the  chill  in  an  ingot  of  9  to  1  reduction  in 
area  (thermal  data.  Appendix  G) . 
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Figure  39:  Composition  versus  distance  from  chill  for  reduced 
section  ingot,  (a)  Theoretical  -  Cq  ■  4.5  per  cent 
copper,  10/1  reduction,  (b)  Calculation  from 
experimental  data  CQ  ■  4.68  per  cent  copper,  9/1 
reduction. 
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the  plates  may  act  similar  to  an  orifice  with  the  maximum  in  velocity 
occurring  slightly  downstream  from  the  reduction. 

The  2/1  reduction  is  the  same  as  in  one  of  the  reduced  section 
ingots  previously  discussed,  Figure  37.  In  fact  the  shapes  of  the 
two  composition  curves  are  very  similar. 

If  comparison  is  made  of  Figures  36,  37,  and  41  it  can  be  seen, 
as  expected  from  Figure  3,  that  the  minimum  composition  is  a  function 
of  the  area  reduction. 

B.  Unidirectional  and  Bi-Directional  Heat  Flow;  Non-Unidirectlonal 
Fluid  Flow 

In  the  theoretical  section  v  was  discussed  as  a  part  of  the 
v  G  X 

x 

parameter,  — —  .  When  v^  ■  0  the  general  macrosegregation  equation 
(20)  reduces  to  that  for  solidification  at  a  chill  face.  This  is  the 
maximum  composition  obtainable  in  unidirectional  solidification  (in 
the  absence  of  melting  and  solidification)  and  for  an  aluminum-4.5 
per  cent  copper  alloy  is  4.92  per  cent  copper.  The  minimum  composi¬ 
tion  is  kCg  ■  .77  per  cent  copper  for  4.5  per  cent  =  which  will 
be  approached  as  v^  approaches  minus  infinity.  For  solidification 

where  both  heat  and  fluid  flow  are  non-directional  v  can  become 

x 

positive  and  local  compositions  greater  than  4.92  per  cent  copper  are 
possible.  In  general,  for  solidification  where  heat  flow  is  unidirec¬ 
tional  and  no  melting  and  resolidif lcation  occurs,  as  v^  increases 
negatively  from  zero,  local  average  composition  will  decrease. 
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1.  Unidirectional  Heat  Flow,  Non-Unidirectional  Fluid  Flow. 

For  unidirectional  solidification  v  refers  to  that  component  of 

X 

fluid  flow  which  is  perpendicular  to  isotherms  (parallel  to  heat 
flow).  In  this  section  three  castings  will  be  considered  in  which 
values  and  the  resulting  compositions  varied  considerably. 

The  first  experimental  ingot  considered  is  shown  in  Figure  42. 

The  "ingot"  comprises  the  bottom  portion  of  the  mold  cavity  (3.8  cm, 
by  15  cm  by  10  cm  high).  The  hot  top  is  as  shown,  feeding  through  a 
neck  at  one  end  of  the  ingot.  Ingot,  neck,  and  hot  top  are  "perfectly" 
insulated  so  all  heat  flow  is  downwards,  out  through  the  water-cooled 
base  plate. 

The  ingot  is  designed  to  simulate  heat  and  fluid  flow  in  solidifi¬ 
cation  of  a  slice  of  a  large  ingot.  Specifically,  the  base  of  the 
ingot  is  intended  to  represent  the  mold  wall,  the  "top"  of  the 
ingot  represents  the  centerline,  and  flow  occurs  left  to  right  from 
the  neck,  as  it  would  downwards  from  the  hot  top  of  a  large  ingot. 

Solidification  in  the  portion  of  the  ingot  near  the  base  is  as  in 
a  simple,  constant  cross  section,  unidirectionally  solidified  ingot. 
This  is  confirmed  by  the  final  segregation  seen  in  this  region  (at 
locations  4.5,  6.1  and  7.6  cm  from  the  chill,  Figure  43).  At  any 
given  location  along  the  casting,  composition  is  essentially  the  same 


at  all  three  distances  from  the  chill. 
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As  the  dendrite  tips  approach  the  Ingot  top  some  differences  In 
flow  behavior  are  to  be  expected: 

(1)  fluid  flow  can  no  longer  be  perpendicular  to  isotherms 
throughout  the  ingot  if  solidification  shrinkage  is  to  be  completely 
fed,  (2)  downward  flow  in  the  vicinity  of  the  neck  must  be  much 
greater  than  that  in  a  classical  unidirectional  ingot  at  a  similar 
distance  from  the  chill.  In  addition,  it  may  be  expected  that  isotherms 
in  the  vicinity  of  the  neck  may  become  crowded  relative  to  the  rest  of 
the  ingot,  and  alter  Qjt\  this,  however,  did  not  significantly  affect 
solute  redistribution  in  the  similar  Ingot  of  Figure  40,  in  which  any 
change  in  movement  of  isotherms  did  not  appreciably  affect  calculate 
segregation. 

Figure  44  shows  the  composition  along  the  ingot  at  locations 
7.9  and  9.8  cm  from  the  chill  (i.e.,  2.1  and  0.2  cm  from  the  casting 
top).  Remarkably  strong  macrosegregation  results, where  compositions 
near  the  ingot  top  vary  systematically  by  as  much  as  .83  per  cent 
which  is  nearly  20  per  cent  of  the  base  alloy  composition  (aluminum- 
4.7  per  cent  copper  alloy). 

Most  notable  feature  of  the  data  is  the  low  solute  region  below 
the  neck.  This  would  be  predicted  from  (2)  above.  In  addition 
the  composition  of  the  ingot  top  at  locations  removed  from  the  neck 
appears  to  be  above  C^;  this  would  be  expected  (in  the  case  of  complete 
feeding)  since  the  flow  here  cannot  be  perependicular  to  isotherms. 

That  is,  vx  is  approaching  zero  and  lies  between  zero  and  the  steady 
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state  value.  This  will  be  discussed  further  In  the  centerline 
segregation  portion. 

The  second  example  considered  in  this  section  is  the  portion  of 
the  previously  described  ingot  that  had  stainless  steel  inserts. 

Figure  45  shows  composition  as  a  function  of  distance  from  the 
chill  for  that  portion  of  the  ingot  shown  which  had  the  plates  to 
restrict  fluid  flow.  The  decrease  in  composition  near  the  chill, 
with  increasing  distance  from  the  chill,  is  as  expected  in  usual 
unidirectional  solidification,  the  minimum  at  2.5  cm,  however,  is 
lower  than  expected;  it  is  not  clear  why  this  is  so. 

It  is  interesting  to  note,  by  comparing  Figures  41  and  45,  that 
the  minimum  in  the  section  with  the  plates  occurs  at  a  different 
distance  from  the  chill  than  minimum  in  the  other  portion  of  the 
ingot.  At  the  location  where  the  right  portion  reaches  a  minimum 
(at  4.40  per  cent  copper)  the  left  portion  is  about  the  same  composi¬ 
tion.  The  left  side  (without  plates),  however,  continues  to  decrease 
to  a  minimum  of  4.16  per  cent  copper-  At  this  location  the  right 
portion  (with  plates)  has  a  local  average  composition  greater  than 
5.0  per  cent  copper.  This  means  that  a  composition  difference  of 
more  than  .8  per  cent  copper  exists  over  a  horizontal  distance  of 
less  than  2.5  cm. 

Near  the  plates  where  is  expected  to  be  0  the  composition  should 
be  the  same  as  the  chill  face  composition  (i.e.,  approximately  5.00  per 
cent  for  the  4.56  per  cent  ■  alloy).  That  fact  that  the  composition 
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just  under  and  between  the  plates  is  higher  could  be  due  to  some 
"exudation"  occurring  as  a  result  of  solid  contraction. 

Just  above  the  plate,  soliditication  also  is  expected  to  occur 
with  vx  ■  0  and  composition  here  is  the  same  as  at  the  chill  face. 
Further  above  this  plate,  a  minimum  is  observed  which  may  be  related 
to  another  set  of  plates  (not  shown)  15.25  cm  from  the  chill. 

The  third  example  considered  is  the  previous  ingot  (9/1  reduction) 

at  one  side  of  the  heavy  section  (Figure  11).  As  discussed 

previously  when  theoretical  calculations  were  made  for  this  ingot, 

velocity  at  a  location  y  (for  a  ^iven  distance  from  the  chill)  was 

assumed  constant  except  exactly  at  the  riser-casting  junction.  In 

the  actual  case,  however,  this  is  not  expected  to  be  true.  As  feed 

metal  enters  the  larger  section  the  absolute  value  of  v  must  be  a 

x 

maximum  at  the  center  and,  just  below  the  section  reduction  is 
expected  to  be  greater  than  the  flow  velocity  in  unidirectional 
solidification,  constant  cross  section;  hence,  final  composition  is 
here  less  than  C^,  Figure  40. 

To  feed  the  heavy  section  out  toward  the  edges,  v^  is  expected 
to  approach  zero  for  x  *  L  where  L  is  distance  from  the  chill  to 
section  reduction.  Hence,  composition  at  the  top  of  the  heavy 
section  (away  from  the  reduced  section)  should  be  that  at  a  chill  where 
v^  is  also  zero  (5.15  per  cent  copper  for  a  4-68  per  cent  *  alloy). 


This  is  seen  in  Figure  40. 
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2.  Centerline  Segregation. 

While  centerline  segregation  is  considered  as  a  separate  problem, 
it  can  be  treated  as  a  limiting  special  case  of  non-unidirectlonal 
heat  and  fluid  flow.  As  discussed  in  the  theoretical  section,  Part 
III,  the  increased  solute  content  at  the  centerline  of  an  ingot  can 
be  related  to  the  magnitude  and  component  of  the  fluid  flow  velocity 
v^.  As  an  example,  consider  the  simulated  centerline  ingot  previously 
discussed,  Figure  42.  In  Figure  46  the  results  of  chemical  analysis 
away  from  chill  are  presented  at  one  location  along  a  direction 
parallel  to  the  heat  flow  direction.  Composition  remains  relatively 
constant  until  the  top  is  approached.  At  locations  close  to  the  chill, 

fluid  flow  to  feed  solidification  is  perpendicular  to  isotherms  and 

8  e 

v  approaches  a  steady  state  value  ■; - -  t  .  As  the  dendrite  tips 

X  1  —  P  Lf 

touch  the  top,  fluid  flow  can  no  longer  be  perpendicular  to  isotherms 
and  vx  approaches  zero.  As  v^  approaches  zero,  local  average  composi¬ 
tion  increases.  The  maximum  occurring  when  v^  =  0  and  composition 
approaches  the  same  as  at  a  chill  face.  That  this  maximum  was  not 
reached  in  this  ingot  at  the  simulated  centerline  was  presumably  due 
to  incomplete  feeding  along  the  very  top  of  the  ingot. 

It  can  now  be  visualized  how  centerline  segregation  can  occur  in  a 
large  ingot  where  solidification  fronts  from  the  ingot  walls  meet  at  a 
central  region.  Before  the  dendrite  tips  touch,  solidification  is  as 
above,  when  v  is  perpendicular  to  isotnerms.  Once  the  tips  touch, 
however,  flow  is  no  longer  perpendicular  to  isotherms  but  can  even 
become  parallel.  When  flow  is  parallel,  *  0  and  local  average 
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composition  will  approach  chill  face  composition  (for  perfect 
feeding).  The  extent  that  local  average  composition  approaches 
chill  face  composition  depends  on: 

1.  How  much  solidification  in  the  volume  element  occurred  with 

v  -  0. 
x 

2.  How  closely  v  approaches  zero  (which  depends  on  details  of 

X 

fluid  flow  as  discussed  in  the  previous  section) . 

As  example  of  this,  consider  the  ingot  solidifying  toward  a 

common  center  from  two  flat  chill  walls,  using  the  experimental  set-up 

described  in  part  I-B,  Appendix  E.  Figure  47  presents  the  results  of 

43 

chemical  analysis  across  this  ingot  at  two  locations  .  Both  show 
increases  in  com*'  ition  at  the  center  even  though  casting  width  was 
only  3.75  cm  an^.  total  solidification  time  was  140  second.  The 
qualitative  similarity  of  the  segregation  curves  in  Figure  47  to  those 
measured  in  large  ingots  is  remarkable;  see,  for  example,  Reference  17. 

Finally,  consideration  can  again  be  given  to  both  the  reduced 
cross  section  ingot  and  the  ingot  with  metal  plates  inserted.  In 
these  ingots,  v^  is  expected  to  approach  zero  at  the  following 
locations,  as  well  as  at  the  chill  face: 

1.  In  the  reduced  section  ingot  at  the  top  of  the  heavy  section 
(see  Figure  40b) . 

2.  In  the  ingot  with  inserted  plates;  on  either  side  of  the 
plates  and  inside  the  plates  (see  Figure  41). 


no 


Local  average  composition  at  these  locations  Increases  in  a  manner 
similar  to  the  top  of  the  simulated  centerline  ingot  and  the  actual 
centerline  of  the  experimental  ingot. 
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IV.  CONCLUSIONS 


1.  Flow  of  solute-rich  liquid  to  feed  solidification  and  thermal 
contractions,  is  the  basic  mechanism  responsible  for  many 
apparently  different  types  of  segregation.  Among  these  are: 

(a)  inverse  segregation 

(b)  centerline  segregation 

(c)  banding 

(d)  negative  cone  of  segregation 

(e)  localized  solute  rich  regions 

(1)  exudation 

(2)  "hot  tearing"  type  segregation 

(3)  solute  enrichment  under  risers  and  hot  tops. 

2.  The  final  composition  at  any  point  x,  y,  z  in  a  casting  or  ingot 
(and  hence  macrosegregation)  can  be  calculated  from  knowledge  of 

(a)  phase  diagram  and  other  physical  constants  of  the  alloy, 

(b)  thermal  history  at  x,  y,  z  during  solidification,  and  (c) 
amount  and  direction  of  fluid  flow  at  x,  y,  z  to  feed  shrinkage 
of  parts  of  the  casting  at  lower  temperatures. 

3.  Laboratory  experiments  on  aluminum-4.5  per  cent  copper  alloy 
qualitatively  and  quantitatively  confirm  analyses  presented. 
These  experiments  comprise  measurement  of  macrosugregation  in 
small  ingots  of  the  following  mold  designs:  (a)  simple 
unidirectional, (b)  unidirectional  with  restricted  fluid  flow, 

(c)  unidirectional  with  variable  heat  extraction  rate  to 
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illustrate  banding,  (d)  unidirectional  with  reduced  cross 
section,  (e)  edge-risered  unidirectional  to  illustrate  centerline 
segregation,  and  (f)  ingot  with  bi-directional  heat  flow  to 
Illustrate  centerline  segregation. 
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V.  SUGGESTIONS  FOR  FUTURE  WORK 

While  much  experimental  work  has  been  done  in  the  area  of 
unidirectional  heat  and  fluid  flow,  further  experimental  work 
seems  desirable  on  non-directional  solidification.  Specifically, 
work  should  be  done  on  "hot  tear"  type  segregation,  segregation  in 
sand  castings,  and  under-riser  segregation. 

In  conjunction  with  this  work  consideration  should  be  given  to 
the  effect  of  porosity  on  macrosegregation  results.  While  a  t  eore- 
tical  and  qualitative  presentation  is  made,  quantitative  resui  ,s 
should  prove  helpful. 

A  new  area  of  investigation  would  be  the  effect  of  magnetic 
fields  on  macrosegregation.  A  sufficiently  strong  magnetic  field 
should  influence  segregation  by  altering  fluid  flow  during  solidifi- 

.4  44 

cation. 

Another  area  needing  study  is  that  concerning  the  details  of 
fluid  flow  in  non-unidirectional  solidification.  Limiting  cases  have 
been  considered  herein,  but  treatment  of  intermediate  areas  (which 
more  closely  approximate  real  conditions)  requires  detailed  fluid 
flow  analysis. 

A  more  applied  area  for  future  work  would  be  application  of  the 
ideas  presented  herein  to  actual  ingot  solidification  processes 
(static,  continuous,  consumable  electrode).  Such  application  would 
require  numerical  analysis  and  clearly  would  require  use  of  high 
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speed  digital  computers.  However ,  the  application  should  be 
feasible,  and  certainly  of  great  engineering  usefulness. 
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APPENDIX  A 


List  of  Symbols 


V  aE 


A*  A’,  Ag,  A^ 


C 

a 


intercepts  of  llquldus  and  eutectic  Isotherm 
movement  equations 


cross  sectional  areas  at  x,  x',  x^ ,  and  chill 
face,  respectively 

local  average  composition  of  liquid  and  solid 
(weight  fraction,  weight  per  cent  solute) 

average  composition  of  the  entire  amount  of 
primary  a  phase  in  the  final  solid  ingot 
(weight  fraction,  weight  per  cent  solute) 

composition  of  eutectic  (weight  fraction, 
weight  per  cent  solute) 

liquid  composition  (weight  fraction,  weight  per 
cent  solute)  within  volume  element 

liquid  composition  (weight  fraction,  weight  per 

cent  solute)  at  x'  when  x_  ■  0 

E 

liquid  composition  (weight  fraction,  weight  per 
cent  solute)  when  feeding  stops 

overall  composition  (weight  fraction,  weight 
per  cent  solute)  in  a  binary  alloy 

average  solid  composition  (weight  fraction, 
weight  per  cent  solute)  after  solidification 


ct 


AC 


81 


8t 


solid  composition  (weight  fraction,  weight  per 
cent  solute)  at  the  liquid-solid  interface 

macrosegregation  (weight  fraction,  weight  per 
cent  solute) 

weight  fraction  liquid 

weight  fraction  solid 

volume  fraction  eutectic 

volume  fraction  liquid 


■ft.  J'.  .  wmm  r- 


'  T"-'« 


~~  'CSJM3 
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8S 

G 

k 

L 


v 


X 


x,  y,  z 


volume  fraction  liquid  at  x'  when  ■  0 

volume  fraction  liquid  when  feeding  stops 
completely 

volume  fraction  solid 


temperature  gradient 

equilibrium  partition  ratio 

total  solidification  length  from  a  chill  or 
wall  (i.e.,  unidirectional  solidification 
L  ■  total  length  of  Ingot,  Ingot  solidification 
from  two  walls  L  ■  width/2  etc) 


slope  of  llquidus  and  eutectic  Isotherm  move¬ 
ment  equation 


(nL  ‘  nE)/nE 


dimensionless  group  of  variables 
time  when  feeding  stops  completely 
time  when  dendrite  tips  reach  L 
non-equilibrium  solidus 
llquidus  temperature 


.  V  (1-8) 

c  6 


velocity  of  eutectic  isotherm 

local  velocity  of  liquid  relative  to  solid 
(resulting  from  solidification  contraction), 
cm/sec 


local  velocity  of  liquid  relative  to  solid  in 
the  x  direction  (perpendicular  to  isotherms) 
resulting  from  solidification  contraction 

coordinate  axles  for  volume  element  under 
consideration 


x',  y',  z '  general  coordinate  axles 

x^  location  of  the  llquidus  Isotherm  when  the 

eutectic  isotherm  is  at  x^,  ■  0 
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location  of  the  llquldus  and  eutectic  Isotherms 
(distance  from  the  chill) 

location  of  "volume  element"  when  feeding  stops 


sol'  ’» 

tion  shrinkage  (p<,  -  PL)/P 

solid;) 

*on  shrinkage  of  eutectic 

rate  c 

perature  change,  #C/sec 

3 

local  average  density  (liquid  and  solid),  gm/cm 

3 

density  of  a  phase  gm/cm 

3 

liquid  density  gm/cm 

3 

solid  density  gm/cm 

densities  of  solid  and  liquid  eutectic 
respectively,  gm/cm-* 
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APPENDIX  B 

* 

Calculations  of  Solute  Redistribution,  Unidirectional 
Solidification  Constant  k.  B 


1.  x  >  x 

c 

Assume  first  the  distribution  of  liquid  in  the  mushy  zone  is  as 
sketched  in  Figure  5b.  Solidification  conditions  are  described  by: 


\  ■  **Lt  + 


XE  "  V  +  aE 


gL  ~  gE  _  1  “  gE 

x’  ’  XE  ’  XL  "  XE 


(Bl) 


(B2) 


(B3) 


Differentiate  (B3)  with  respect  to  time  at  constant  x' 


d  t 


1  -  gF  dx  x'  -  x  . 

■  - ( »  d r  +  < >  Tt  (xl  -  *e» 


(B4) 


Rewrite  (B3)  for  x'  ■  x,  g£  *  gL;  differentiate  with  respect  to 
time  at  constant  g^.  Set  equal  to  -e/G  for  linear  g^  and  T  versus  x: 


<  —  > 
1  dt  ’ 


8L  "  8E 


3x, 


i  -  gc  Kfr  )(xl  -  V  -TT 


£ 

G 


(B5) 


Now,  for  constant  p^»  equation  (22)  is  written: 


*LVX 


dx. 


B  r  r  “°L  .  ,  E  . 

1  -  B  ^  dt  X  'gE  dt  ' 


(B6) 
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Substitute  (BA)  In  (B6)  and  Integrate: 


x-x_  dx_  .  x-x_  .  dx_ 

gLVx  "  "  (  T-6  )1(1_8E)(xt  -x  ^  dt ~  +  1  (x ,  -Xp^dt  (xL~xE^  +  8E  dt"  J 

L  E  L  E 


(B7) 


Combine  (B5)  and  (B7) : 


v  G 
x 


x-Xp  dXp  .  x-Xp  dXp 

,  1  (1-gE>  1  dT  +  I  (l^>  dt  (VXE>  1  +  gE  dT 

1-B  <8,)  g,-8r  -  3x. 


,  >L  eE  „  3  , 

(  - )  «r—  (x  -X_)  +  T7" 

l-g_  3t  L  E  3t 


....  (B8) 


Substitute  In 


x  - 


X 


L 


bV 


b-1 


(B9) 


(BIO) 


dt  (XL  -  V  'If  <]tL  ■  XE)  ■  b  (nL  •  nE)tb"1  <B11) 


n 


n.  -  np 
l  E 

n. 


q  “E 

and  (B8)  becomes 


(B12; 


v  G 
x 


1  ( 

2  '  1  -  6 


B _  )(  gL(gLna  +  2b)  +  "a8E' 


gL<gL"q  +  b) 


(B13) 


where  b  -  1  -  g£  -  nqgE 

Solute  redistribution,  and  from  this,  macrosegregation  can  now  be 
determined  analytically  by  substituting  (B13)  in  equation  (20)  and 
integrating  from  ■  1  at  CL  ■  CQ  to  gL,  CL<  The  analytic  expression 
so  obtained  is  cumbersome.  For  small  g^,  however,  (B13)  reduces  to: 
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gTnn  +  2 

)  [  “k-9 -  ] 

*L"q+  1 


(B14) 


which  is  also  the  expression  obtained  by  assuming  the  distribution  of 
liquid  in  the  mushy  zone  as  given  in  Figure  5c.  Substitution  of 
typical  values  for  unldirectionally  solidified  aluminum-4.5  per  cent 
copper  shows  that  results  obtained  from  (B14)  are  at  most  a  few  per 
cent  different  from  those  obtained  using  (B13) ,  and  (B14)  is  therefore 
used  for  subsequent  analysis. 


As  example  of  the  difference,  consider  the  case  of  k  ■  .172, 

v  G 

0  ■  .055,  n  ■  .52.  Calculated  values  of  obtained  from  the  two 

q  e 

equations  are  given  in  Table  Bl: 


TABLE  Bl 


v  G 

x 

Values  of  -7—  for  Aluminum-4.5  Per  Cent  Copper 


v  G 

_(E_Q .  B 1 4)_ 

v  G 

-7-  (Eq .  B 1 3)  _ 

Per  Cent 
Difference 

£  1 

.0489 

0.0480 

-  1.6 

.0522 

.0517 

-  1.0 

.0567 

.0578 

+  1.9 

.0569 

.0582 

+  2.3 

.0902 

(eutectic) 
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Now,  substituting  equation  (B14)  in  equation  (20),  and  integrating 
from  the  limits  gL  ■  1  at  cL  *  CQ  to  gL»  CL»  the  final  solute 
redistribution  expression  is: 


C 


L 


-(1-k) 


C0gL 


giV  + 1 

n  a  +  1 

q 


»-«>  »-fc) 

a 


(B15) 


where  a 


1 


6 

2 


2.  x  <  x 

c 


For  calculations  of  solute  redistribution  at  locations  close  to 
the  chill  face  a  procedure  similar  to  the  foregoing  is  followed  except 
that  for  x  <  x  and  g.  >  gt,  g.v  (see  equation  B6  above)  is  calculated 

C  L»  Li  L  X 

by  integration  of  equation  (22)  from  the  limits  0  to  x,  not  xr  to  x. 

£ 

The  resulting  expression  for  g  *  0  is: 

Li 


V  X  ,  8  ■  2(11  -  V  +  gLno  (x  -  2xE> 

e  "  2  '  1  -  8  '  ,  ,2,  . 

(x  -  xE)  (gLnq  +  1) 

where 

X  *  gLaE("o  -  1}  +  Vl 

XE  p  n  +  1 

L  q 

Substitute  (B17)  in  (B16) 


(B16) 


(B17) 


v  G  0 

x  x  .  B 

e  2  k  1  -  B 


2p  +  (gT Q  -  x)n 
_ t _ 9. 

gE-2 


(B18) 
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where 

P  *  nq*  -  Vq  +  ®E  '  aL 
Q  *  "qX  2Vq  +  2aE  ’  2aL 


Substitute  (B16)  in  equation  (20) 


1  dCL 


d8. 


where 


1  -  k  CL  N  +  M*l 


6x(2P  -  xn  ) 

N  -  - 3- 

2 

2P 


M 


xBQn  +  2(1  -  6)P' 
2P2 


(B19) 

(B20) 


(B21) 


(B22) 


(B23) 


Integrating  between  the  limits  gL  -  1  at  ■  CQ  and  gL, 


C 


L 


N  +  MgL 
N  +  M 


1  -  k 
M 


) 


(B24) 
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APPENDIX  C 

Calculation  of  Solute  Redistribution  in  Final  Transient 
Region  During  Unidirectional  Solidification 


From  Figure  (C-l) 

tv 

d  -  xy)  -  -  /  vxdt 

t° 

V 

for  n  *  0,  x.  -  x_  -  constant 
(]  L  L 

e  dxg 

Vx  ■  1  -  $  dt 


v 

where  ■  fraction  liquid  at  xv 
Differentiating  (C3)  and  rearranging 


(Cl) 


(C2) 


(C3) 


dXv  ■  <XL  '  XE)  d« 


V 

L 


(C4) 


Substituting  in  (Cl) 

t 

1  -  XV  "  f  -  TT"ii  Idxv  '  (xL  ‘  XE)d(!L  1  (C5> 

t° 

V 

V  V 

at  time  ty,  xy  -  xy,  gL  -  gL 
at  time  t^,  xy  -  L,  g *  - 


1 
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Substituting  these  limits  in  (C5) 


L  -  x 


8 

1-6 


v 

/  dx 


~  t  -  v* *1 

1 


Integrating  between  limits 


1  -  x 


v  "  l  -  g  Kxv  "  L>  “  “  1)J 


Solving  for  g£  yield* 


v  1  V 

6r  ■  1  -  T. - T~ 


<±  -  V 

(\  -  *,3)6 
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4-J 

a  a> 

T>  O  ^ 

•h  y  cx 


O  SL 

ca  u 


E 

o 
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APPENDIX  D 

Calculation  of  Solute  Redistribution*  Unidirectional 
Solidification;  Constant  k,  6;  Pyramid  Shape 

Ingot 


1.  x  >  x 


From  equation  (23)  when  g£  ■  0 


AgLVx  "  "  T“-~B 


dgL 

— -  x'A' 
dt 


(Dl) 


where 


A*  -  (1  -  £  )  A0 


(D2) 


from  (BA)  with  n  ■  0,  x.  ■  x_  ■  constant 

q  L  E 


dgi  1  dxE 

— “  -  (  - - -  )  f  — -  1 

dt  1  ^  -  xE  ’  1  dt  J 


(D3) 


Substitution  of  (D2)  and  (D3)  in  (Dl) 


(1  -  L)gLVx 


X 

/ 


dx. 


1  -  B  '  xT  -  x_  dt 
L  F. 


(1  - 


r)Sx' 


(D4) 


integration  yields 


6  -  <  1 


1  -  B  v  L  -  x  7  dt 


dx  (x  -  x  ) 
)  TT  [L  ,r  E 


(xL-  XE) 


2  2 
x  - 

2(xL  "  V 


(D5) 


from  (B9)  when  gg  ■  0 


x  -  x. 


XL  -  XE  L 


(D6) 
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From  (B17)  when  n  ■  0 

q 


XE  *  X  +  aTgL 


(D7) 


where  aT  ■  aL  - 


(D8) 


From  (B5) 


_e 

G 


(D9) 


Substitution  of  (D7)  and  (D9)  in  (D5) 


2L(xl 


-  XE)  +  2xaT  *  4gL 
2(xL  ‘  XE) 


) 


2.  x  <  x 

c 

Integration  limits  on  (Dl)  change 

X  dg' 

-  •  i4-e  !  i r  3x'A'  <D10) 

0 

Substitution  of  (D2)  and  (D3) ,  integration  and  substitution  of 
(D9)  yields 

1*L  _  _  _ B_  ,  _ L_  \  r  x(2L  -  x)  , 

e  1-B'L-xm  2gL  1 


(Dll) 
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APPENDIX  E 

Experimental  Casting  Procedure 


A.  Introduction 

1.  Unidirectional  Heat  Flow. 

The  typical  experimental  set  up  shown  In  Figure  E-l  was  used, 
except  for  one  series  of  centerline  segregation  experiments,  for  all 
castings  and  Ingots  made  as  part  of  the  current  Investigation.  Since 
unidirectional  solidification  is  desired,  essentially  all  heat  must 
be  extracted  from  one  surface  of  the  ingot.  This  was  accomplished  by 
doing  all  casting  and  solidification  In  an  electrically  heated 
furnace. 

After  placing  the  entire  assembly  in  the  furnace  with  coolant  pipes 
extending  outside,  the  unit  was  heated  to  600-700#C.  When  the  plaster 
reached  the  desired  temperature  (usually  1-2  hours  after  the  furnace 
temperature)  the  mold  was  filled  through  a  funnel  extending  from  the 
outside  of  the  furnace  into  the  cavity.  The  assembly,  with  liquid 
was  held  for  several  minutes  to  allow  convection  to  subside.  The 
coolant,  either  air  or  water  (depending  on  the  desired  solidification 
rate)  was  turned  on  and  the  ingot  solidified.  With  the  four  sides 
and  top  insulated,  heat  removal  for  solidification  was  through  the 
bottom  surface,  resulting  in  fully  columnar  structures. 

Casting  and  Ingot  designs  were  varied,  as  discussed  in  Experimental 
Results  and  Discussion  section,  to  produce  desired  macrosegregation 
effects.  However,  for  casting  solidified  with  unidirectional  heat  flow 
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the  top  was  always  open  to  the  furnace  atmosphere  and  temperature, 
the  bottom  surface  rested  ov  a  chill  through  which  the  heat  was 
extracted,  and  all  other  surfaces  were  insulated  with  plaster. 

2*  Centerline  Segregation  Experiment. 

The  centerline  segregation  experiment  was  designed  to  produce 
two  unldirectionally  solidifying  fronts  meeting  at  the  ingot  center, 
Figure  E-2.  Casting  procedure  was  the  same  as  in  the  straight  uni¬ 
directional  ingot  with  water  or  air  used  as  coolant.  (Piping  is  not 
shown  in  schematic  drawing,  however,  continuous  flow  was  supplied  to 
both  chills.) 

B.  Mold  Preparation 

In  many  ingots,  thermocouples  were  used  to  obtain  continuous 
records  of  solidification.  In  all  cases  these  were  chromel-alumel 
couples  with  output  continuously  monetered  on  a  L  &  N  Speedomatic  W, 
Multipoint  (12)  recorder. 

When  It  was  necessary  to  locate  many  thermocouples  accurately, 
the  pattern  for  the  ingot  was  first  shaped  from  foamed  polystyrene. 

The  thermocouples  (alumina  sheathed,  tips  coated  with  boron  nitride), 
were  placed  in  the  styrofoam  pattern  at  the  desired  locations.  The 
assembly  was  placed  inside  a  stainless  steel  flask,  the  flask  was 
then  filled  with  liquid  plaster  which  was  allowed  to  set.  The  result 
was  that  the  thermocouples  remained  accurately  located  since  they  were 
held  in  place  by  the  plaster. 
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The  mold  was  next  fastened  to  the  chill  and  placed  In  the  furnace 
with  the  foamed  polystyrene  still  in  place.  Upon  subsequent  heating 
the  polystyrene  vaporized  and  left  a  deposit-free  cavity. 

lien  castings  were  made  without  thermocouples,  wood  patterns  were 
used  rather  than  styrofoam.  The  greased  wood  pattern  was  placed 
inside  a  flask,  liquid  plaster  poured  around  the  pattern  and  allowed 
to  set.  After  setting,  the  pattern  was  removed  and  the  mold  fastened 
to  the  chill.  The  assembly  could  then  be  placed  in  the  furnace. 

C.  Metal  Melting  Procedure 

All  melts  were  prepared  from  aluminum-4.5  per  cent  master  heats. 
These  had  been  previously  cast  into  5  to  7  pound  ingots  from  high 
purity  aluminum  (99.99  per  cent)  and  50-50  aluminum-copper  OFHC  master 
alloy.  All  melting  was  done  in  "Fiberfrax"  coated  graphite  crucibles 
using  a  gas-fired  furnace.  After  preheating  the  crucible  for  an  hour 
or  more,  the  charge  was  melted  and  superheated  to  a  temperature  of 
690-720°C.  After  stirring  the  melt,  a  chlorine  degas  was  carried  on 
for  10  minutes.  After  stirring  again  for  one  minute  a  gas  sample, 
under  reduced  pressure,  was  taken.  When  the  sample  was  gas-free, 
the  crucible  was  removed  from  the  furnace,  a  chemical  analysis 
specimen  taken  up  into  a  7  mm  vycor  tubing,  and  skimmed.  Pouring 
temperature  into  the  mold  was  about  680  to  700°C.  After  pouring, 
another  chemical  analysis  specimen  was  taken  from  the  remaining  liquid 


metal. 
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Figure  E-l:  (a)  Schematic  diagram  of  unidirectional  solidification, 

mold-chill  assembly. 

(b)  Cross  section  of  assembly. 
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Figure  E-2:  (a)  Schematic  diagram  of  centerline  segregation, 

mold-chill  assembly. 

(b)  Cross  section  of  assembly. 


APPENDIX  F 


Chemical  Anlysis  by  Fluorescence 


A.  Introduction 

When  any  material  Is  bombarded  with  electrons  of  X-rays  of  high 
enough  energy  a  characteristic  line  spectrum  is  emitted.  This  can 
be  used  as  a  basis  for  either  qualitative  or  quantitative  chemical 
analysis.  For  qualitative  analysis  the  various  elements  in  a  sample 
are  made  to  emit  their  "characteristic  lines",  these  elements  may  then 
be  identified  by  analyzing  the  emitted  radiation  and  showing  that  the 
specific  wavelengths  characteristic  of  the  element  are  present.  By 
comparing  the  intensity  of  these  lines  with  the  intensities  of  lines 
from  suitable  standards,  a  method  of  quantitative  analysis  is 
available.  For  multicomponent  systems,  comperison  with  standards  is 
complicated  by  interactions  between  elements.  When  dealing  with 
binary  systems  however,  it  is  a  rapid,  accurate  method  for  obtaining 
reproducible  chemical  analysis. 

In  the  present  work  a  General  Electric  XRD-5,  X-ray  Spectrometer 
was  used  with  a  tungsten  target  tube.  The  primary  white  radiation 
from  the  tube  was  used  to  fluoresce  aluminum-copper  samples  on  an 
area  1.25  cm  by  1.8  cm.  This  sample  size  insured  that  the  analysis 
was  representative  and  no  microsegregation  effects  are  recorded. 

That  is,  we  are  measuring  local  composition  by  integrating  over  many 
dendrites,  which  are  typically  lO-lOOOp  center  to  center  spacing. 
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The  secondary  radiation  from  the  aluminum-copper  alloy  is  then 
diffracted  from  a  single  crystal  of  known  d  spacing  (lithium  fluoride, 
d  spacing  2.014) . 

The  Bragg  law  predicts  that  for  each  angular  setting  of  the 
crystal, radiation  of  only  a  single  wavelength  is  reflected.  The 
spectrometer  was  set  at  45.08°  which  for  the  crystal  and  spacing  used 
reflects  the  copper,  wavelength.  The  intensity  of  this  "character¬ 
istic  line"  was  measured  and  compared  to  a  standard  intensity  curve. 

B.  Standards  Preparation 

40 

Standards  were  prepared  with  the  same  apparatus  used  by  Bower  in 
a  study  of  solidification  in  thin  plate  castings.  (Shown  schematically 
in  Figure  F-l) .  By  drawing  a  vacuum  and  placing  the  molds  in  the  melt 
molten  metal  is  forced  between  the  two  copper  chills,  spaced  1/8  inch 
apart,  and  allowed  to  solidify.  Compositional  variations  were  obtained 
by  successively  adding  high  copper  master  alloy  to  low  alloy  heats. 

A  10  pound  heat  of  alumlnum-1  per  cent  copper  was  melted  in  a  gas 
furnace  using  a  graphite  crucible.  After  degasing  with  chlorine,  samples 
were  sucked  into  the  chill  mold  at  1500°F.  Master  alloy  of  50  per  cent 
copper-aluminum  composition  was  added  to  increase  alloy  content  of  the 
melt  to  about  2.5  per  cent  copper.  Again  the  melt  was  degassed  and 
samples  sucked  into  the  chill  mold.  By  progressively  adding  alloy,  a 
series  of  standards  were  prepared.  After  cooling  to  room  temperature 
the  samples  were  sanded  on  coarse  paper  to  remove  any  exuded  layer. 

They  were  then  sanded  down  to  300  grit  emergy,  a  finish  smooth  enough 
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for  use  In  fluorescent  analysis.  Table  F-l  lists  the  compositions 
obtained. 


TABLE  F-l 


Standard 


Composition-%  Copper 


A 

B 

C 

D 

E 


1.12 

2.69 

4.35 

5.83 

8.48 


After  sanding  the  cast  surfaces  on  the  edge, two  pieces  weighing 
about  one  gram  each  were  removed  from  opposite  sides  and  used  for 
wet  analysis.  The  average  of  these  compositions  was  taken  to  be 
the  composition  of  the  plate. 

C.  Standard  Curve  Determination 

Fluorescence  of  the  five  standards  was  carried  out  as  previously 
described.  The  intensity  from  Standard  E  (maximum  copper)  was  measured 
first  with  the  power  adjusted  to  give  approximately  100,000  counts  per 
second.  This  approaches  the  limits  of  the  scaler  used  for  counting. 

At  greater  than  100,000  counts/sec  the  unit  can  be  flooded  and 
inaccuracies  result.  After  the  intensity  from  Standard  E  was  determined, 
Standards  A  through  D  were  measured.  Since  the  intensity  of  copper,  Ka 
is  being  used  for  comparison,  decreasing  copper  content  means  decreasing 
intensity.  After  all  five  samples  are  measured  Standard  E  is  normalized 


TT 
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to  I/I  ■  1.  The  Intensity  of  standards  is  A  to  D  are  divided  by  I 
(the  intensity  of  Standard  G)  to  give  normalized  Intensities.  (These 
values  are  plotted  in  Figure  F-2) . 

D.  Composition  Determination 

When  the  Intensity  from  a  series  of  unknowns  is  measured,  Standard 
E  is  also  measured.  The  intensity  of  the  unknown  divided  by  the 
intensity  of  Standard  E  (Iq)  gives  a  fraction  less  than  one  (if  the 
composition  of  the  unknown  is  less  than  8.48  per  cent  copper).  This 
fraction  can  be  directly  translated  to  a  per  cent  copper  from  Figure 
F-2. 


Standard  practice  was  to  check  one  of  the  Standards  A  to  D  each 
time  a  series  of  unknowns  was  determined  to  be  sure  reproducibility  of 
results  was  obtained.  As  a  further  check,  periodically  samples  which 
had  been  analyzed  by  fluorescence  were  also  checked  by  wet  chemical 
methods.  Agreement  was  found  generally  to  be  within  +  0.05  copper. 
Another  method  of  insuring  accuracy  was  the  redetermination  of 
standard  curves.  Three  different  times  during  the  investigation  the 
standard  curves  were  redetermined.  Each  time,  no  significant 


difference  was  found  between  curves. 
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To  Aspirator 


Figure  F-l:  Schematic  diagram  of  apparatus  used  to  prepare 
standard  chemical  analysis  specimens. 


Figure  F-2:  Standard  chemical  analysis  curve. 
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APPENDIX  G 

Outline  for  Composition  Calculations 

The  following  is  an  outline  of  the  method  used  to  calculate 
chemical  composition  at  various  locations  in  unidirectionally 
solidified  ingots,  using  thermal  data  and  relationships  developed 
in  the  theoretical  section.  As  example,  the  calculations  are  given 
for  a  point  7.9  cm  from  the  chill,  in  the  reduced  section  ingot 
described  in  the  Experimental  Results  and  Discussion  section. 

A.  Thermal  Data 

Temperature  versus  time  curves  at  the  different  thermocouple 
locations  are  made  directly  from  the  output  recorded  on  the  L  &  N 
Speedomatic  W  Recorder  (Figure  G-l) . 

B.  Temperature  Versus  Fraction  Liquid 

A  temperature  versus  fraction  liquid  curve  is  drawn  assuming 
g_  ■  0  and  fraction  liquid  distribution  in  the  mushy  zone  is  linear 
(Figure  G-2a).  For  comparison,  a  curve  of  actual  fraction  liquid  in 
the  mushy  zone  is  shown,  Figure  G-2b  (see  Macrosegregation  Versus 
Distance  from  Chill  in  Theoretical  Section). 

C.  Calculation  of  v 
- x— 

By  using  the  data  from  Figures  G-l  and  G-2a,  Figure  G-3  is 
constructed  to  be  used  for  graphical  integration  to  find  vx  as  follows: 
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(1)  At  a  constant  x,  Ag  Ax  Is  determined  by  adding  up  the  number  of 

L 

squares  from  g.  ■  0  to  g_  ■  1  between  t  and  t  .  -  and  multiplying  by 
l  l  n  n+1 

the  area  of  each  square.  Area  to  the  left  of  the  Ingot  reduction  Is 
weighted  9  times  that  to  the  right  (ratio  of  area  reduction) .  That  Is 
In  equation  form. 


(2)  Equation  (23)  Is  solved  for  the  case  of  no  eutectic: 


A28LVx 


1-81 


/ 


x„ 

£. 


K 


t  . .  -  t 
n+1  n 


dx'  + 


K 


1  -  8  "2 


Cn+1  ~  tn 


dx* 


(Gl) 


Equation  (G-l)  Is  written  for  the  reduced  cross  section  Ingot, 


for  x  >  x  and  for  x_  >  0  where: 
c  E 


c 

X 


distance  from  chill  to  reduction  In  cross  section 

distance  from  chill  to  location  where  composition  is  being 
calculated 


Ag^  ■  change  in  fraction  liquid  at  x  from  t^  to  tR+^ 

A^  ■  area  of  lower  (large)  section 
A2  *  area  of  upper  section 

For  the  general  case  of  constant  cross  section  at  locations  in 
the  ingot  unaffected  by  the  reduction,  equation  (G-l)  becomes: 


gLVx 


1-8 


( 


AgLAx 


t  . ,  -  t 
n+1  n 


) 


(G2) 
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D.  v  Versus  b. 

-x —  ■  nL— 

Figure  G-4  is  obtained  by  taking  as  being  linear  and  averaged 
over  the  interval  AgL>  That  is  the  determined  value  is  plotted  at: 


8i 


*8L^n  +  ^8L*n+l 


(G3) 


where  g^  ■  fraction  liquid  at  x  at  time  t 


t  +  t  . . 
n  n+1 


This  is  plotted  as  Figure  G-4. 


E.  Calculation  of  G/e 

G/e  is  calculated  by  first  replotting  the  data  of  Figure  G-3  as 
in  Figure  G-5  and  then  the  reciprocal  of  these  slopes  are  G/e  at  a 
given  x  for  various  g^.  In  Figure  G-5  sample  curves  for  two  values 
of  g^  are  shown.  In  practice  a  whole  family  of  curves  is  necessary 
to  calculate  G/e  as  a  function  of  g^. 

F.  G/e  Versus  b^ 

Figure  G-6  shows  G/e  as  a  function  of  g .  Using  this  and  Figure 
v  G  L 

X 

4.  —  as  a  function  of  g.  can  be  calculated. 

*  e 

G.  Calculation  of  Local  Averaae  Composition 

v  G 

From  the  steady  state  solute  redistribution  expression  —  -  is 

found  as  a  function  of  C.  .  C.  being  calculated  from  the  values  of  gT 

L  L  _  T. 

v  G 

x 

assuming  steady  state  values.  The  final  plot  is  — versus  lnC^ 
which  is  necessary  to  solve  the  solute  redistribution  (equation  20). 
This  is  shown  in  Figure  G-7. 
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Average  composition  Is  found  as  follows: 


Solve  by  taking  two  cases  of  macrosegregation  equation: 


(1  -  k) 


% 

8t 


v  G  dC 

-  (i  -  b>u  +  -?-> 


(G4) 


also  take  case  where  v  -  0,  at  chill  face  where  maximum  segregation 

X 

can  be  calculated 


(1  -  k) 


dg, 

8t 


dC. 


-  (1  -  6) 


(G5) 


Subtracting  (la)  from  (1) 

g,  CL  v  G  dC 

(1  -  k)  In  — 7  ■  (1  -  B)  /  -r-  7-*  (G6) 

Bj  c  Lt 

co 

v  G 

.  A  plot  of  — —  versus  InCj  as  in  Figure  G-7  allows  graphical 
integration  to  solve  the  general  macrosegregation  equation. 

The  area  under  the  curve  will  give  the  term  inside  the  integral  of 
equation  (3): 


area 
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Figure  G-4:  Calculated  vx  versus  g^  for  the  9  to  1  reduced 
section  ingot  at  x  ■  7.9  cm  from  the  chill. 


Figure  G-l:  Thermal  data  for  9  to  1  reduced  section  ingot  (composition  versus  distance 

from  chill  presented  in  Figure  38).  Numbers  are  location  of  thermocouple 
(cm  from  chill). 
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Figure  G-5:  Constant  fraction  liquid  curves  plotted  as  a 
function  of  time  versus  distance  from  the 
chill. 
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Figure  G-6:  G/e  versus  fraction  liquid  (g^)  at  7.9  cm  from  chill 
in  the  9  to  1  reduced  section  ingot. 


